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PREFACE 


111 tile present work our nienn efforts have been to 
study the static properties of light nuclei (A.k.40) which are 
obtained fron Hartree-Eock calculations with realistic inter- 
actions. I'liis i7ork mainly consists of five chapters. Chapter I 
is mainly an introduction giving a brief account of the indep- 
endent particle model of the nucleus, and the Hartree-Pock and 
Brueclcaer-FIa.rtree-Eock approximations of solving many-body 
problem. Several Hartree-Pock and 3rueckner-Hartree-Pock 
calculations that have been done in the past are also summarized 
and the problems taken up for the present vfork have been 
mentioned. Chapter II describes the nucleon-nucleon inter- 
actions which can be used in nuclear structure studies. ?fe 
have chosen the following two interactions for our purpose; 
effective interaction of Shakin, ¥aghmare and Hull derived 
from Yale potential by using a unitary -model-operator approach 
and the Sussex interaction of Elliott et. al. obtained directly 
from experimental phase shift da.ta by using an auxiliary potent- 
ial method. Chapter III contains the Hartree-Pock theory in 
detail and the matrix method of solving the Hartree-Pock 
equations. Ihe formulae used in the calculations of various 
ground state properties are also given. Two important cor- 
rections due to the coulomb energy*' and centre-of-mass motion 

of the nucleus are estimated. Both these corrections are 
included in a self-consistent manner in our calculations, 

U) 



unlike tiie esirlier Hartree-Fock calculations of various other 
authors. Ihese corrections are estimated by those authors 
externally after solving the Hartree-Fock equations. Single 
oscillator approximation limit to the Hartree-Fock solutions 
has also been discussed. In Chapter IV we present the result 
of spherical Hartree-Fock calculations for the nuclei ^He, ^Be 
^^0 and ^*^Ca. The properties calculated are the binding 
energies, the single particle energies, r.m.s. radius and 
charge radius, matter and charge densities, spin-orbit split- 
tings , etc. It is found that the binding energy results, 
obtained with effective Yale interaction for all the nuclei 
considered are in satisfactory agreement with the correspond- 
ing experimental values if all the second-order terms of the 
Yale potential are accounted for. Sussex matrix elements givd" 

better results for the binding energies in the first order 

40 

compared to the first order in Yale potential. For Oa, 
however, Sussex matrix elements give binding energy in excess 
of the experimental value. Rest of the properties calculated 
are in fair agreement v/ith the experimental data. The results 
obtained in the single oscillator configuration approximation 
to Hartree-Fock wave functions are also presented. Single 
particle wave functions and self-consistent single particle 
potentials obtained in the Hartree-Fock calculations with 
effective Yale interaction have been studied in detail. The 
Hartree-Fock single particle wave functions for ^^0 have been 
compared with the conventional harmonic oscillator wave 
functions. The general nature of the wave functions is the 



same for all the nuclei. The Hartroe-Pock single particle 
potential is non-local and state— dependent . The static limit 
to this non-local potential has "been discussed. 

Pinally , in Chapter Y wo have studied the energy 
18 18 10 

levels of Oj P and ^Be in the Hartree-Pock representation. 
It is assumod that the nucleons move independently of each 
other in a Hartree-Pock self-consistent field and the low- 
lying states of these nuclei may arise due to the excitations 
of only the last two neutrons to the higher states. The single 
particle energies and wave functions describing the behaviour 
of these tv/o nucleons are obtained from Hartree-Pock calcul- 
ations. The nucleon-nucleon potentials used are the effective 
Yale interaction of Siiakin et. al. and Sussex interaction. It 
is seen that the HP spectra are compressed compared to experim- 
ental levels in all the cases. Probable reasons for this 
discrepancy have been discussed. 



CHAPTER I 


INTROPUCTIOIT 

I-l Independent Particle Model 

IHe study of the nucleus is a many-hody prohlem which 
cannot he solved exactly. Instead, studies of complex nuclei 
have been carried out by various phenomenological models. There 
are several nuclear models which can be broadly divided into 
two classes. First class consists of so-called ‘strong inter- 
action models’ which treats the nucleus as an assemblage of 
closely coupled particles, e.g., liquid drop model, Bohr and 
Mottelson model, etc. These models have been successful in 
explaining many collective properties of nuclei but we shall 
not discuss them here. The second class is of Independent 
particle models'’’ which assumes that the nticleons move independ- 
ently of each other in an average potential. Recent evidence 
suggests that the independent particle model provides a better 
description for low-energy phenomena though the coupling of 
nucleons by virtue of their mutual interaction can by no means 
be neglected. The most successful of these independent particle 
models are the shell model, the optical , model and the Nilsson 
model. The uniform potential is spherical for shell model, 
complex for optical model and spheroidal for Nilsson model. 

Shell model, proposed by Mayer in 1948 and developed since 
then, has been very successful in explaining some regularities 
in the binding energies of nuclei, spin and parities of some 
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of the low-lying nuolear levels, oertain multipole transitions, 
P-deoay, etc* The optical model was primarily designed to 
explain the nuclear reaction phenomena and it successfully 
accounts for the scattering cross-section data of low and medium 
energy reactions. Nilsson model can predict well the ground 
state spins, moments of inertia, energies of excited state ro- 
tational hands, equilibrium deformations, and vibrational pro- 
perties of deformed nuclei, which are not described well on the 
basis of simple shell model. In all these independent particle 
models the nucleons fill up the successive lowest quantum levels 
of the mean potential well and move in definite orbits. Ihis 
picture provides the qualitative features of many nuolear pheno- 
mena. To get the quantitative agreement with experiments the 
residual interaction between the particles is introduced. It 
mixes up the orbits to some extent but the basic character of 
the model is essentially the independent particle motion. 

With increasing addition in our knowledge of nucleon- 
nucleon force from the study of two, three and four nucleons 
systems, it was thougiit that there should be a way of correlat- 
ing the average potential field of the nucleus with- the nucleon- 
nucleon interaction. The strong objection to this was the fol- 
lowing. Most of the realistic nucleon-nucleon pptentials 
which fit the two-body data contain a repulsive core. Such a 
repulsive core will give rise to strong short-range oorrelatio.tis, 
which are not described well by an independent particle model. 
However, there are many experimental evidences in favour of 



independent particle notion, e.g. the nuclear reaction experi- 
ments, in which an incident particle knocks out a nucleon from 
target nuclei, show that single particle structure does prevail 
for the nucleons inside the nucleus. Besides, elastic scatter- 
ing shovB re sonance s which are explained in terms of optical 
potential model. Ihus, the fact that independent particle model 
is successful in explaining a large number of experimental data 
implies that short-range correlations do not play an important 
role in many nuclear properties and the nucleons inside the 
nucleus behave- to a certain extent as independent particles in 
a common potential. Recent advances in the theory of nuclear 
matter by Moszkowski and Scott (Moi 60) and G-omes et, al. 

(Gos 58) have explained the mechanisms which reduce the effect 
of strong interaction on the nuclear motion. In particular, the 
essential features of Moszkowski- Soott method could be described 
as follov/s: 

(i) The scattering of two nucleons inside a nuclear medium 
is characteristically different from the collision of free 
nucleons. The reason for this is that inside the nuclear 
medium Pauli principle forbids scattering to states which are 
already occupied and thus, the number of available final states 
is limited. The nucleons far apart can interact through low 
Pourier components of the force. Such interactions are gene- 
rally forbidden inside the nucleus since the momentum transfer 
is insufficient to raise the particles above the Perm sea. 

The short-range part of the potential corresponds to. high 
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momentum transfer and scatters the particles to unfilled states. 
Thus, beyond a certain healing distance the relatiTe wave 
function of two nucleons approaches that of noninter acting 
particles. This is in sharp contrast to scattering of free 
nucleons where the wave function never heals up and is phase 
shifted relative to free particle wave function even at infi- 
nity. The healing distance is of the order of IT* and the ave- 
rage separation of the nucleons in nuclear matter is approxi- 
mately 1.66P. Consequently, the effect of two-body correlations 
would be small. Besides, the fact that healing distance is 
smaller than the average separation implies that the probability I 
that more than two particles will simultaneously be at distances | 
smaller than the healing distance would be small and three-body 
correlations would be unimportant. I 

(ii) The effect of the attractive part of the potential on 
the nuclear phase shifts is opposite to that of the repulsive 
part. Thu^ a part of the attractive potential can be combined^ 
with the repulsive core so that the two together give zero phase : 
shifts. Consequently, -(lie phase shifts due to the strong nucleon- 
nucleon interaction are the same as thaf due to weak attractive I 
long-range part of the force beyond the healing or cut-off [ 

distance. 

As a consequence of these two effects particles move ! 
almost independently with very weak interactions. This implies 



that the nuclear Hamiltonian H of a system of A nucleons can 
be written as 


H = + V, 


( 1 - 1 ) 


where is the independent particle Hamiltonian and v is the 
residual interaction. I'he average potential field, which each 
nucleon feels, must be some sort of average of the interaotions 
it has with all the other nucleons. The residual interaction 
between two nucleons must be the interaction left over after 
this average effect is removed. A systematic way of determining 
Hq, starting from fundamental nucleon-nucleon interaction, is 
the self ‘-consistent Hartree-Hock (HE) method. The ground state 
of the nucleus is represented, in the first approximation, by 

■T 

a nondegenerate eigenstate of H^ i.e. 




( 1 . 2 ) 


where is the lowest eigenvalue of H^. We shall briefly 
describe the HE method of determining^^ and E^. 

^ 2 Hartre e-Eoc k Me thod 


The Hartree-Eock (EE) method consists in assuming tha,t | 
the wave function^^ of the nucleus is an anti symmetrized pro-/ 
duct of single particle wave functions. To determine the best | 
single particle wave functions there are many equivalent critexia 
Eor example, one can minimize the total energy or one can require 
self-consistency between the single particle potential field of 
II and the average potential seen by the particle as a result 
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of its interao tion \=o.tli ail the other nucleons. The single 
particle wave functions which minimize the total energy are the 
solutions of following HI equations 

<ct|T||3) + <a|Tllfj> = , (1.3) 

where a and p specify the single particle states, T is the 
kinetic energy and U the HP single particle potential which is 
related to the nucleon-nucleon potential 7 hy the following 
equation, 

■ A 

<afujp^ = ^<aYtv|PY>» (1-^) 

Y=1 

Thus, the HP single particle potential U is calculated in terms 
of Y hy a sum over occupied single particle states, but the lattej 
cannot be obtained without a prior loiowledge of U. This is the 
HP self-consistency problem and is solved by an iteration 
interpolation procedure, starting from a trial guess and improv- 
ing the wave functions until self-consistency is achieved. The 
HP field may be spherical or deformed depending upon the nuclei 
studied. Thus, for example, the HP field is spherical for the 
doubly closed shell nuclei "^He, ^^0, ^^Oa and ^^^Pb. The ground ; 
statei^Q is a Slater detemtinant wherein all the single particle 
states are occupied up to Permi level e^,, Por . deformed nuclei 
single particle- states are described in a deformed HP field. 

The ground stateQ^ is again a Slater determinant in which all 
the deformed orbitals are occupied up to Perm! level. 
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The HF method not only desorihes the independent parti- 
cle picture of the nucleus hut can also he used in the discus- 
sion of collective motion of the niicleons, e.g.,the vihration 
and rotation of the nucleons as a vi/hole can. he explained hy let- 
ting the self-consistent HF field become time— dependent. Besides, 
the pairing correlations of the two-hody force can he tahen into 
account in an extended W theory, Icnown as Hartree-Foch- 
Bogoliuhov-Valatin theory. This method extends the variational 
principle to a class of trial wave functions of independent 
excitations of quasi particles which include pairing correla- 
tions. In this case the total Hamiltonian H can he written as 

H = + v*^^ , (1.5) 

OF 

where is the independent quasi particle Hamiltonian and 
OF 

v“ is the residual interaction between these quasi particles. 

HF theory has been successfully applied in atomic and 
solid state physics. In fact, the idea of self-consistent field 
was first intutively introducod hy Hartree and was successfully 
applied to electronic structure of atoms, later this was 
extended hy Foch who included the exchange effects and estab- 
lished the HF equations using a variational method. Thus, the 
HF method in nuclear physics is essentially borrowed from 
atomic physics, hut the applica.tion of this method to nuclear 
structure calculations is much more tedious than in the case 
of atoms because of the following reasons; 
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(i) The nuol eon-nucleon interaction Y is not known 
uni quely. 

(ii) Most of the realistic nucleon-nucleon potentials which 
fit the ti/vo-hody data contain a hard core. The matrix element 
of such a potential in a basis of independent particle wave 
functions is infinite, 

(iii) Complete HE calculations are very time-consuming and 
cannot be carried out mthout high-speed computors. 

Inspite of all these difficulties several spherical 
and deformed HI' calculations have been done during the last 
few years. A detailed survey of all these calculations is not 
possible here, Por an excellent review of deformed HP calcu- 
lations we refer to an article by Eipka (Eia 68). Since our 
main interest lies in the spherical HH calculations, we shall 
briefly review only such calculations that have been carried 
out in the past, Sphorical HE caLoulations by Davies et. al. 
(Das 66) have been ca-rried out using simple velocity- dependent 
force acting only in relative s-state. They have studied the 
nuclei ^^0, ^^Si, "^^Ca and ®®Zr and obtained qualitative 

agreement for varioj^s cclbulated ground state properties with 
experiment. In the calculations of Muthukrishnan et. al. 

(Mun 65 ), Yamaguchi potential has been used and the results 
obtained for various nuclei are very similar to those of Davies 
et. al. (Das 66). Abgrall and Monsonego (Abl 66) have used a 
central force with a nonsaturating exchange mixture, j\1s®, 
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tiiere are calculations done by Krieger et. al. (Krr 66), whose 
effective interaction includes three-body contact term propor- 
tional to cube of density. In their later calculations, Davies 
et. al. (Ian 68) have employed a more complete two-body poten- 
tial including both even and odd noncentral potentials and . 
giving a fair fit to two-body data. Ihey have studied in detail 
a number of closed shell nuclei in the range 16 4-1^208. 

Although, in some of these calculations the calculated HF proper- 
ties are in reasonable agreement with experiments, -the major 
drawback of these calculations is that they were carried out 
with nucleon-nucleon potentials especially designed for use in 
HF calculations, Ihe parajneters of these potentials were • 
adjusted to fit the nuclear matter data so that they gave a 
fair fit for the binding energies of various nuclei (particu- 
larly heavy nuclei) with the experimental data, but they did 
not fit the two-body data. In fact those effective potentials 
do not contain any features of the realistic nucleon- nucleon 
potentials, such as hard core or nonlocality, one-pi on- exchange 
behaviour at large distances, etc. First time a realistic poten- 
tial (labalcin potential) in the HF calculations has been used 
by Kerman and c-O -workers (Ken 66). These authors have studied 
the nuclei ^^0 and in considerable detail and find that 

labakin potential fails to give adequate binding in first-order 
and gives too large a value for the spin-orbit splittings. 

Later on, Pal et. al. (Pal 66) have improved the agreement for 
the binding energy by calculating the contribution of second 
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order term. In all tile HE* calculations described above the 
nucleon-nucleon interaction ?/as smooth and did not contain any 
hard core singularity. The first theoretical foundation for 
the use of singular potential in HP calculations has been laid 
by Shakin et. al. (Shn 66 , S7a)- Iheir theoiry is based on the 
unit ary-model- operator approach of Villars (Vis 63) and is 
applicable to any general singular interaction. Ihe fact that 
short-range correlations do not play an important role in many 
nuclear properties implies that independent particle HP wave 
functions are significant and can be obtained from an effective 
interaction from vfhich the effect of hard core has been removed. 
This was the philosophy adoiited by Shahin et. al. in dedx|icing 
a smooth effective interaction for use in HP calculations. 

Thus, all the problems from (i) to (iii) have now largely been 
solved. There is another approximation method proposed oy 
Brueckner and oo~workers for dealing with the problem of hard 
cores- We describe it very briefly below. 

1.3 Brueckner Theory 


In the Brueckner theory motion of a pair of nucleons 
is described by means' of a correlated wave funotionJ^pC^,^) 
which vanishes inside the core region | ^-r^ | ^ r^ . and asympt- 
otically approaches the unperturbed wave funotion(^^p: 

( 1 . 6 ) 
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'Jis tile solution of a modified form of two-nuoleon Soiirodinger 
equation known as Brueckner-Bethe-Goldstone equation (Bee 57), 
and is related to () through a reaction matrix (K-matrix) defined 
as 


• U.7) 

Because of the strong short-range correlations, generated hy 
the two-nucleon potential, it is assumed that the motion of the 
pair will not he affected hy the other nucleons through any 
multipartiole clustering, the other nucleons only modify the 
energy of the interacting pair through HF potential, and Pauli 
principle forbids scattering to states which are occupied hy 
other nucleons. This is Iniomi as independent pair approximation. 
The interaction between the nucleons inside the nuclear medium 
is determined hy a reaction-matrix which differs from the two- 
hody potential over the short-range of correlation distances in 
the nuclear wave functions. The defining equation for K-matrix 
i s 


K = V + V f K » 


( 1 . 8 ) 


Q is the Pauli operator which takes care of exclusion principle 
and e is the energy denominator which takes into account the 
binding effects of average nuclear field. In a basis of indep- 
endent particle wave functions cpjs together with associated 
eigenvalues EJs, the matrix elements of E can, he written as 
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K 

a 


p , y6 , CTT , y6 


0T 


(1.9) 


The sura over a and x is to le carried out over empty states only. 
The single particle energies octsuaarliig in equation (1.9) are • 


E = 
a 


(^) 

^ Zm^aa 


+ 


> (K, 

y( occ ) 


— K ") 

ay, ay ay,ya‘' 


( 1 . 10 ) 


TEe second term in equation (1,10) defines tlie single particle 
potential U, i.e. 

U = / (K -K ) . (1.11) 

aa ^ ay, ay ay,ya^ / 

y(ooc) 


In the limit of a weak interaction K can be replaced by the 

potential V and the single particle potential defined by . i 

equation (1.11) is the HE potential. Starting from a given \; 

set of eigenstates 9^ (^ means n^^ iterate) the new set 9^'*'^ I 

obtained from the single particle potential will in general ^ 

not agree with the input. The problem of obtaining agreement 

is' the self-consistency problem of HE. In the Brueckner theory , 

a new self-consistency problem, arises since a change in the 

n n+l 

representation 9^->9^ not onl.y changes the single particle 
potential, but also through the change in the energy spectrum 
changes the reaction-matrix. This in turn reacts back on the 
single particle potential and eigenfunctions. This is the 
Brueckner self-consi stency problem of obtaining self-consistent 
energy spectrum and reaction-matrix. In infinite nuclear 



matter the eigenstates are plane waves and only Brueokner self- 
consistency problem exists. In finite nuclei Brueckner and HF 
self-consistency problems are coupled together and exact 
Brueckner-Hartree-Fock (BHF) calculations become quite compl- 
icated. Besides, in finite nuclei exact treatment of Pauli 
operator causes some difficulties and various approximation 
methods are used to calculate the reaction-matrix. Many authors 
have bypassed the HF self-consi stenoy problem by using pure 
■ oscillator wave functions (Hog 67, 68; Kor 66, 6?; Mar 66). -Such 
calculations encounter only Brueokner self-consi stenoy problem 
and are called Brueckner- type calculations. In his calculations 
Wong (Wog 67) has studied qualitatively the ground state prop- 
erties of ^*^0a, ^^0 and ^He using a number of hard core potent- 
ials. It is found that realistic potentials underbind these 
nuclei by 2-4 MeV per particle, a result similar to that 
obtained in pure HF calculations with realistic potentials 
(Ken 66: Shn 67b), Some BIIF calculations for finite nuclei 
have also been reported in the literature. First few BHF 
o.alculations (Brr 58, 61; Man 63: Kor 65) have been done in 
coordinate representation, but none of them are exact. Calcul- 
ations by Brueckner and coworkers have been done in the local 
density approximation. In tiiis approximation, the reaction- 
matrix elements for finite nuclei are taken to be simple 
averages over matrix elements for nuclear matter, calculated for 
the density at the centre- of-ma.ss position of the two inter- 
acting nucleons. In the calculations carried out with Gammel- 
Thaler potential, light nuclei are found to be underbound by 
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2-4 MeV per particle but discrepancy decreases with, increasing 
A. However, the equilibrium nuclear radius becomes too small 
for heavy nuclei. The results for ^^Ca, ^^Zr and ^*^®Pb 

are summarized in references (3rr 61) and (Man 63)* BHF calcul- 
ations in an oscillator basis have been recently reported by 
Davies et. al» (Das 69). As in the pure HD calculations, the 
single particle wave functions are expanded in terms of harmonic 
oscillator wave functions and the calculations have been carried 

out with Hamada- Johns ton potential for the nuclei ^^0, '^^Ca, 

4-8 2 08 

Ca and Pb. Iheir main conclusions are that nuclei are too 
small and Underbound. 

¥e thus see that HP is very much alive even when the 
nucleon-nucleon interaction is shrgT-ilar, if not the simple HP 
something more sophisticated like BHP is valid- In fact, it has 
been shown (Grg 68) that the effective interaction of Shakin 
et. al., obtained by the method of canonical transformations, 
can also be expressed in terms of K-matrix. However, in the HP 
calculations c.arried out with this effective interaction 
( Shn 67b) there is no Brueckner-self-consistency involved since 
in the calcula,tion of second-order terms in the effective inter- 
action these authors have used an average value for the energy 
denominator e and angle averaged Pauli operator Q, Por the 
purpose of present work we come back to simple HP calculations 
which do not involve any Bruckner- self-consistency problem. 



Pro'blems for Present Worlr 

The following problems have been undertaken in the 
present work; 

(i) Spherical in' calculations with effective Yale inter- 
action have been c.arried out for “Sle, ^Be, ^^0 and ^^0. This 
interaction has been widely used in several spherical and 
deformed HF calculations of various nuclei. Spherical HP 
calculations for ^^0 and "^^Ca with this interaction have been 
done by Shakin et, al. ( Shn 67b). Present work is an extension 
of their calculations to other lighter nuclei. ¥e have also 
included coulomb and centrc-of-mass corrections in our calcul- 
ations in a self-consistent -mannex whereas in the calcul- 

ations of Shakin et. al. the effect of these corrections on 
the total ground state energy was estimated externally. The 
effect of these corrections on various ground state properties 
ha's also been studied in detail. The results are presented 
in Section. IV. 2 of Chapter IV, 

(ii) Elliott et. al, at Sussex (Elt 68a) have recently 
p ublishe d relative matrix elements of the nucleon-nucleon inter- 
action in an oscillator basis obtained directly from experim- 
ental phase shift data. As a test of their matrix elements 
Elliott et. al. (Elt 67) have aalculated the spin-orbit split- 
tings of light nuclei (A$41) and obtained reasonable values. 

A number of deuteron (Elt 68b) and triton and alpha particle 
properties (Jan 69) have also been calculated and satisfactory 



agreement mth experiments lias been obtained. The success of 
those calculations suggesis that Sussex matrix elements are 
relia-ble ^nd can be used in the calculations of more complicated 
nuclei. We have used these matrix elements in the spherical 
calculations for the nucle-i "^He, ^^0, ^^0, and The 

results are given in Section IV, 3 of Chapter IV, A detailed 
comparison of the Sussex matrix elements results with those 
of several other oalouln-tions, particulanly with those of poten- 
tials obtained from a fit to experimental phase shifts, has been 
made in Section IV, 4 of Chapter IV. Since the HP formalism is 
well-known to give a good account of the ground state properties, 
our calculations also serve as a test of the saturation prop-, 
erties of Sussex matrix elements. 

(iii) We have calculated the low-lying excited states of 
10 1 R 1 R 

■^^Be, 0 and P using HP v/rave functions. As in the pure shell 

model calculations, it is assumed that low-lying states of these 
nuclei arise due to the interaction of last two nucleons and 
remaining nucleons are treated a.s inert. The single particle 
parameters needed in the calculation are obtained from self- 
consistent HP calculations. 

Several shell model spectroscopy calculations for the 
energy spectra, of these nuclei have been reported in the liter- 
ature. In most of the eanlier calculations the residual inter- 
action bet^jyeen the last two nucleons v;as taken to be purely 
phenomenological, e.g. some simple exchange force with G-aussian 
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or Yukav/a radial dependence. Ihe adjustable parameters of the 
potential were determined from best fit to experimental data on 
one or more nuclei. Such calculations have been done by Waghmare 
( ¥ae 64 ) s , Elliott and Elov^ers (Elt 57)»and many others, 

The residual interaction so obtained depends upon the config- 
urations taken into account and the energy levels which are 
included in the calculation and cannot be compared with any 
standard realistic potentials. A realistic interaatien in the 
energy level calculations for A = 18 nuclei was first used by 
Kuo and Brown (Kuo 66). These authors? starting from Hamada— 
Johnston potentials have obtained an effective interaction which 
when used in the standard spectroscopy calculation, gives a good 
description of the energy levels of these nuclei. Olement and 
Baranger (Git 68) using TabaMn potential have obtained results 
for A = 18 nuclei very close to those of Kuo and Brown. In all 
tho calculations described above pure harmonic oscillator wave 
functions and experimental single particle energies have invari- 
ably been used. The use of harmonic oscillator wave functions . 
is preferred because in this representation separation between 
relative and centre-of-mass coordinates can be made exactly and 
two-body matrix elements of the nucleon-nucleon potential can be 
easily 'calculated. Very recently, Kahana (Kaa-Preprint) has rep- 
orted energy level calculations for A = 18 nuclei using Woods-Saxon 
vrave functions. Ko attempt rill now has been made to obtain the 
single particle energies and single particle wave functions 
from some self-consistent calculations. HE calculation gives 



US a. picture very close to siiell model, i . e., it provides us with 
a self-consistent single particle level spectrum and associated 
with each level is the single paxticle wave function. These 
single poxtiole parameters serve as a "basis for any shell model 
spectroscopy calculation and the residual interaction would "be 
same as that used in the self-consistent HF calculation to 
obtain the single particle parameters. We have carried out the 
calculations with effective Yale interaction and Sussex inter- 
action. The details of the ca,loulations and the results are 
presented in Chapter V. 


In. a somewhat different approach followed by some 
authors one tries to dispense with any potential picture and the 
two-body matrix elements which enter any shell model spectroscopy 
calculations are treated as arbitrary parameters and adjusted to 
give best possible fit to spectra of several nuclei. Such 
calculations have been done by Arima (Ara 68) for oxygen isotopes 
^^0, ^^0. In his calculations ^^0 is regarded as inert 
core a,nd low-lying states of these nuclei are assumed to arise 
from neutrons in (Id^y'g configuration. X" ~ ^ matrix 

elements have been calculated from least square fit to these 


levels. ^ = 0 ma,trix elements have been calculated by extend- 
ing the calculations to cover the nuclei F, ^F and he as 
well. Similar calculations have also been done by Federman 
and Tolmi (Fen 65). It would be interesting to compare the 
empirically determined two-body matrix elements with those 
calculated from a realistic potential for this configuration 
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space. Sucli a comparison has heen made in Chapter Y, 

In Chapter II vre describe various types of nucleon- 
nucleon interactions available for use in any nuclear structure 
calculations. The prescription used to obtain the effective 
Yale interaction from Yale potential by Shakin et. al. and 
Sussex matrix elements directly from phase shift data by Elliott 
et. al. has also been described briefly. Tha details of the HE 
theory and various formulae used in the calculations are given in 
Chapter III. The results of HE calculation are presented in 
Chapter IV. The results obtained for the energy levels of 
various nuclei obtained v/ith HE v/ave functions and different 
types of residual interactions are given in Chapter V. 

All calculations have been carried out on IBM 7044 
computer at Indian Institute of Technology, Kanpur. 



CHiPTEE II 


HUCLEOl-IUGLEOE lETERAGTIOES 

Tile exact nature of the nuoleon-nuoleon Interaction is 
still not known and it is not possible to present it in a- close 
deductive form. However, the nucleon-nucleon interaction in 
its nonrel ati vi sti o version can be represented in terms of a 
potential. The nucleon-nucleon potentials, available for use 
in nuclear structure calculations, can be broadly classified _ 
into following categories: 

w 

(i) Realistic potentials 

(ii) Effective potentials 

(iii) Potential matrix elements obtained directly from 
phase shifts. 

¥e shall describe these potentials in the following sections. 

11*1 Realistic Potentials 

Realistic potentials are obtained directly from an 
analysis of the deuteron and two-body scattering data. The 
available data (Marr 68: Eos 67) oonsistof in pp and np scatter- 
ing: total cross-sections, angular distributions, and different 
types of polarization measurements and also s-wave scattering 
lengths. The two-body bound state data are deuteron binding 
energy, quadrupole moment, D-state probability, etc. In the 
energy range 0 to 300 MeY there are about one thousand pieces 
of data left to be fitted, after discarding quite a few of them 
considered imrealistio. There are two ways of handling the 



problem. One can either tiy to fit the data as they oome from 
the experiments or one can use a phase shift parameterization of 
the data and compare these phase shifts with those extracted 
from potentials. Ihe latter procedure is generally preferred. 
Phase shift analysis of the experimental data has been carried 
out mainly by tv/o groups; one at Yale led by Breit (Brt 62) and 
the other at Livermore (Marr 68). Signell and co-worhers 
( Sil 65) have also made a number of useful analyses. Only one 
set of phase shifts emerges and the phases are fairly well 
determined. It is to these phase shifts that the potential 
parameters are fitted in most of the models. No attempt is made 
to fit the data at energies above about 300 MeV. At higher . 
energies pion production and other relativistic effects become 
important and two-nucleon Schrbdinger equation is inadequate. 
However, the energies of nucleons in collision in nuclei are 
expected to be of the order of 160 MeV or so. Thus, we have 
reasonably accurate infoimation about most of the partial waves 
which would play a dominant role in nuclear phenomena. 

The long-range part of the interaction is rather 
unambiguously Icnown from meson- theory . It is determined by the 
one-pion exchange potential (Hun 57) corresponding to the light- 
est meson which can be exchanged between two nucleons. The pion- 
nucleon coupling constant is measured from pion experiments or 
by nucleon-nucleon scattering data. However, at the pion rai^e 
other terms corresponding to the exchange of heavier mesons, 
such as p , , 'n-mesons (Bm 67; Lon 68), dominate one-pion 



exoliange potential (OPEP) by a factor of two or tiiree. Again, 
the one boson exchange strength (OBEP) is deterisiined by the 
nucleon-nncleon scattering or alternatively by calculating the 
me son- nude on coupling constants, for example from SU^ symmetry. 

The nature of the short-range part of the interaction 
is not exactly known. The most general two-body interaction is 
a non-local one, which, to be translational invariant, must have 
the form 

= 6(1-?/) V(rjf"0 . (2.1) 

The relative and oentre-of-mass coordinates are defined by 

'r = i^-r^, ^ = • 2 (^+r‘^) etc. (2.2) 

Meson-theory indicates that the correct two-body interaction 
may be non-local at short-range, but most of the potential 
models are local potentials, for which 

V(r}r^O = 6(^r’/) V(rO . (2.5) 

II. la' Local Potentials 

The first group of local potentials contains a repulsive 
core of infinite magnitude — the so-called hard core. The 
assimption that nucleon-nucleon potential becomes repulsive at 
short distances was first made by Jastrow (Jaw 51). It was 
strongly supported by the two-nucleon scattering data through 
the change in sign of the ^SQ-state phase shift from its low- 
energy positive sign (attractive potential) to a negative sign 



23 


oharaoteristio of a repulsive potential above 200 MeV. [Dv/o suoh. 
most widely known potentials are of Yale group of Breit and col- 
laborators (Laa 62) and of Hamada and Johnston (Haa 62). Yale 
potential is of the following form; 


(2-4) 

where x is the internucleon distance measured in units of pion 
Oompton wave length p,i.e. 


X = p ^r 


(2.5) 


( 2 ) 

Y^ ‘'(x) = one-pion exchange potential 


= + ^ 1-^2 6 (=^)} . ( 2 - 6 ) 


- ,, 2'J X 
X 


12 ^ 


where ^ and ^ are the spin and isospin vectors for the two 
nucleons s ^2 is the usual tensor-operator defined as. 


A 

s. 


(^.,1) (^..^) - (^.-7p) , 


12 , 2^^1 


r- 2- 


(2.7) 


^ — fo- = 0-94- 
411 2^^2 14 _ 


for singlet even states, 


1_ 

14 


otherwise. 


( 2 . 8 ) 


Here (i=G» 1? t S, q) are the central, tensor, spin-orbit 
and quadratic spin-orbit potentials defined by 



'^jj_(x) ‘= +cO for , 

\(x) = for X >x^, (2-9) 

x^ is 0.35. Tile ra<3i'a3 of the hard core is 

. r^ = O. 5 II 6 F for singlet even and triplet odd* 

= 0,5002F for singlet odd and triplet even 

states. (2.10) . 

The parameters of the potential can he found in reference 
(Mar 66). In some cases the idealized hard core is replaced hy 
softer repulsive core which is quite strong hut not infinitely 
large. The examples of soft core potentials are those of 
Bressel, Kerman and lomon (Brl 65) and of Reid (Red 68). Outside 
the repulsive core BKL potential has nearly the same form as that of 
Hamada-Johnston (HJ ) with a few changes in parameters. 

The next group of local potentials is typified hy the 
boundary- condition model of Reshhaoh, lomon and Tuhis (leh 64). 

In this the one-plon and two-pion exchange potentials are used 
outside a radius r^ at which an energy dependent boundary condi- 
tion is specified to fit the high-energy scattering data. 

Another group consists of meson- theoretic potentials 
which t ah e into account exchange of several bosons ^ in . . 
addition to m-mesor.^ such as p »U}s Ypand o-mesons. Different 
combinations of masses and coupling constants give a variety 



of potentials. Recent attempts in developing such potentials 
have been done hy Biyan and Scott (Bm 67) and G-reen and col- 
1 ah orators (G-m 67)* 

II, Ih Non-local Potentials 

A non-local potential is a non-diagonal operator, which 
implies that the' potential energy is a matrix in coordinate 
space. The Schrodinger equation then has the form, 

- ^ + jvCr^r’O ^(^0 dP = E , (2.11) 

where mis the nucleon mass. Because of hermitioity, V('?‘}^') 
must he a symmetric matrix, real in the energy range O^^OO 
MeV. Such potentials are developed hy Yamaguchi (Yai 54), 
Tahahin (Tan 64), Mitra and co-workers and others. Yamaguchi 
potential is a simple separa.hle potential of the form, 


V(^l^) = - Au(r) U(^0. (2.12) 

Tahakin potential is a series of separable terms of Yamaguchi 
type and is of the form; 


V(r!?0 


\ y ^-rC-s(3^)s(^0+h(r)h(r/)]'(||gj(r)'l^^Tgj(^0 


m. 


ith lK2 


(2.13) 


2 

Here ^ vector spherical harmonios. 


defined in terms of the spherical harmonics y^(^) and spin 
states A.g by 





(2.14) 


(tSm^,mgiJM) is the Glebsh-G-ordbn coefficient defined on page 37 
of reference (Eds 57). In expression (2.13) hecause the sum- 
mation is 'done over t and it allows E-mixing, i.e., especially 
S-D mixing. ' 2 .^ is the isospin projection operator. Ihe g-term 
is attractive and h-term is repulsive. Iaba3ri.n has given the 
form of g and h in momentum space* This potential gives only 
a qualitative fit to the two-body data. Moreover, it does not 
have the one-pion exchange feature for the long-range part of 
the potential’. Thus, the labaMn potential is not very real- 
istic but it has been very useful in nuclear structure oalcul- 
-ations because of the ease with which its matrix elements can 
be calculated, 

II-,lc Velocity-Dependent Potentials 

Such potentials are of the form 

V = V(r, p) + constraints for hermiticity, (2.15) 

where p is- the momentum operator. The most oommonly used 
velocity- dependent potential is that due to Green (Grp 62) and 


is of the form 


W(r) and U(r) Have the Gaussian radial dependence. 

II. 2. Effective Potentials 

• It. is well Imovm that the approximation in 'whicH one 

considers nuclear particles moving in some average potential 
field and .interacting with a v/eak residual interaction. . . is . 
Highly successful in tHo' s'Jrudy of nuclear properties.. ..In an 
independent particle model calculation one usually . considers ■ 
the effects of only few shells, IHe neglected shells renorm- 
alize the interaction in the space of. .the shells considered in 
one’s calculation* Eenoinalization of the nucleon-nucleon 
interaction also takes. place due to admixture of veiy high 
..energy, orhitals into the nuclear wave functions, This results 
from strongly ..repulsivp nature of the nuclear force at short 
distances . and also from specific charaoterLstic of tenser 
force *• In fact, the singular nature of the nucleon-nucleon 
interaction at short distances gives rise to strong short-range 
correlations .which- are not properly descrihed hy an independent 
particle, .wave function. Thus, in the nuclear structure calcul- 
ations one must. use the unmodified or the free nucleon-nucleon 
interaction in, a configuration space described by wave functions [ 
which properly take into account short-range correlation effects.} 
..iltematively, one can use a modified or effective interaction, j 
•in -a truncated- configuration space consisting of ..independent | 
pa,rticle wave functions,* ■’ Thus, the problem -is to find an 
effective interaction which would give the same results 

when used with 'unoorrelated wave functions in a limited _ f 





configuration space, as the true interaction V would give -nth 
correla,ted wave fane ti ons in an unliinited configuration space, 
I’or many years, in some simple minded calculations was 

chosen in an empirical manner and only a vague idea existed 
about its possible dependence on the choice of sin^e particle 
orbitals and configuration space and its relationship to the 
realistic ti^ro-body interaction. It was only after the develop- 
ment of Brueckner theory (Brr 5 5) that several different appro- 
aches to the relationship of and V were given. In the 

Brueckner theory the nucleon-nucleon interaction Y is replaced 
by an effective interaction, commonly known as K-matrix and is 
defined in equation (1.8). Several different approaches to 
relationship of Y^^^ and V, based on Brueckner theory, have 
been given by vanious authors, lo name a few are the separation 
method of Moszkowski and Scott (Moi 60), reference- spectrum 


method of Bethe, Brandow and Petschek (Bee 63) 
method (Edn 59), etc. 


Eden-Emery 


The reaction-matrix differs from the two-body potent- 
ial over the short-range of correlation distance in the nuclear 
wave functions. There is an alternative way of introducing 
the effect of short-raage correlations in the effective inter- 
action, namely the method of canonical transformation. This 
approach (which in essence is close to Moszkowski- Scott 
technique) has been followed by Shakin, Waghmare and co-workers 
(Shn 67a). We shall briefly describe below their method of 
obtaining the effective interaction. In their method* the 



correlations in the wave function are introduced via a unitary- 
model-operator such that if the long-range properties of the wave 
function, are specified by a. function ^ ^ ... the 

corresponding correlated state is , 

^ ... ^2 ••• * (2.17) 

An effective Hamiltonian is defined in the space of the uncor- 
related functions as 


H 


eff 


e , 


(2.18) 


where H is the true nuclear Hamiltonian. Since a transformation 
of the form described in eq_uation (2.18) leaves the eigenvalue 
spectrum unchanged, H and H^^^ arc completely equivalent with 
respect to energy cad-culations. But,, besides the desired two- 
body termSjHg^^ contains many-body terms induced by the unitary 
transformation. Carrying out the expansion of H^^^ and neglect- 
ing three and many-body terms, one finds that the effective 
Hamiltonian is, 


H 


eff 




n^ng 


n^n2n^n^ 


X' Wh’ l“3V 


3 

(2.19 ) 


where H^^^ has been expressed in the second-quantized notation. 
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The operators and a are the usual oreation and annihil- 
ation operators for the single particle states. The suh- 
soripts n^ and refer to the quantities (n, h , m) 
necessary to specify the orbitals for particle motion in 
harmonic oscillator potential, T is the kinetic energy 
operator atid \2 is the nucleon- nucleon interaction which 
may or may not ha,ve a hard core. The potentials and U 2 
are used to generate the single particle wave functions. 

The form of U is gcneraJ-ly taken to be of harmonic oscillator 

1 2 

type to simplify the calculations, i.e., U = kr . The 
operator S is assumed to introduce only central short-range 
two-body correlations in which case equation (2,17) becomes 




( 2 . 20 ) 


where^’^s are the uncorrelated wave functions and(p’s are the 
oorr^ted wave functions occurring in equation (2.19). 7^2 

can be written as a sum of t\TO pa.rts, v ^2 which is diagonal 
in relative orbital angular momentima I and v^^ which is off- 
diagonal in Ij, i.e. 


= V + 

V12 + V^, 


( 2 . 21 ) 


Clearly v^"^ receives contributions only from tensor part 

of the nucleon-nucleon interaction. v ^2 is separated into 

p e 

a long-range part v ^^2 s. siiort-range part v ^2 such that 

the short-range part produces no energy shift in the pair 
state. This is a natural extension to finite nuclei of the 
zero .phase shift condition of Moszkowski and Scott (Moi 60) 



for nuclear matter calculations. Ihe main effect of v| 


OD 

T 


is to admix higli-momentum orbitals into nuclear wave function. 
This admixture is not readily trea.ted by the separation method 
used for v^^* Thusj the tensor correlations are treated on 
a different footing than the correlations due to short-range 

S \ 

part of separation of into 

achieved by solving the following two equations s 


( T^+T 2+Ui+U 2 +^ 12 ^^ n^n^ 


and (2i+'^2^""'l‘*'^2'^yn^n2 


niHg 

( 2.22) 

= (e +£ )(1'; , 

- Un iio^^nTn^ 

( 2 . 23 ) 


£ and £ are the energies of single particle states n, 

II 2^ XI 2 

and respectively. For hard core potentials such as the 
Yale or the EJ potential it is necessary to introduce a short- 
range pseudopotential (VP) to carry out this separation 
procedure. In that case 


^2 = ^'^12 + + ^^12 ■ 


(2.24) 


and one requires that 




"r"2 


. (2.25) 

1 -2 - ""1""2 


In general? therefore, one has 

... 


H 


eff < '^i t ® 1 ^ 2 '} ®h. 


^1^2 


1 ^2 


^1^2^3^4 




<$n,n, - y If n,.,) 


(2. 26 ) 



If tile shell-model or Hartree-fock oalculations are limited 
to orhitals having only fairly low-momentum components, almost 
all of the contribution of v^ to the binding energy and to 
the effective interaction would be missed. Most of the 
contribution of v^^ may be tnicen into account by renormalizing 
the effective Hamiltonian in the independent pair approxii^ 
mation (Kuo 65) s 


H 


eff 




^ 1^2 






<fn^n2 I (^*-2 - W) h- ^“^0“ | 


If the second- order term in the pseudopotenti a,l 
included, a term VP tP should be added to 
(2.27). 


(2.27) 


is al so 
in equation 


The prescription described above is quite general 
and is applicable to any singular interaction. Using 
expression (2.27), Shakin et. al. (Shn 67a,67b) have evaluated 
the matrix elements of the Yale potential for various states 
of relative angular momentum. An attractive square-well 
type of pseudopotential (VP) has been used for those states 
of relative motion where the interaction is repulsive and 
also in Sq and states (to keep the separation distance 

const.ant at 1.01). lor the treatment of intermediate states 
in the calculation of second-order terms the approximation 
of Kuo and Brown (Kuo 65) is used, which consists of using 



plane-wave intermediate states and an angle averaged operator 
to take aoGOunt of the Pauli principle in the intermediate 
states. A simple parameterization of the energy denominator 
e is made- Ihe details of these caloulations can he found in 
references (Shn 67a, 67h), Results for the relative matrix 
elements of effective Yale interaction are tabulated (Shn 67b) 
and will be used in the present work. 


It can be, seen the. t the unitary-mo del -ope rat or 
method of obtaining the effective interaction is in some 
sense analogous to Brueckner theory. In fact it is possible 
to find a unitary transformation (Grg 68) for which the 
transformed two-body potential is 

Ygff = |(K + K'^) (2.2S) 


Expansion of defined in equation (2.18) is 


Heff = , (2.29) 


^(n) 


where ^ is aJi n-hody operator* In first-order 


Heff = H 


(1) jj(2) ^ 


(?.30) 


when many-body terms are neglected. Applying the Hartree- 
Eock formalism to in equation (2.30), one finds the 

same expression for the ground state energies as in the 
first-order of Brueckner theory, except for the replacement' 
of K by the hermitian combination ^(K+K'*’), The higher-order 
corrections in both cases anise due to many particle effects, 



more aoourately, the moJiy-tody terms in the result of a 
unitary transformation have their counter par is in the cluster 
corrections of Brueckner theory. 

There are quite a few other effective interactions 
obtained directly from realistic intera.ctions, for example, 
the effective interaction of Ktio and Brown (Kuo 66) derived 
from HJ potential and later improved by Kuo (Kuo 67). Their 
effective interaction includes the contribution of renormt‘.,.. 
alization effects due to core-polarization and has been 
widely used in nuclear- structure studies. To name a few 
other effective interactions used in the literature, are 
those of Wong (V/og 68) obtained from HJ potential and of 
Clement and Baranger (Git 68) derived from Tab akin potential. 
These effective interactions have been obtained following 
the lines of Brueckner theory. 

Of late, there has been a tendency to develop 
potentials primarily for tho use in nuclear structure calGul“ 
ations. The principal aim of these works is to determine 
an effective potential having no hard core singularity and 
giving the correct energy, density and perhaps the symmetry 
energy in nuclear matter. Only a moderate fit to the two- 
body data is required, in additional requirement is that 
second and higher-order contributions to the energy of 
nuclear matter be small in order to improve convergence of 
perturbation expansion. These potentials are effective 
potentials in some sense and it is presumed that they are 



not very diiferent from the realistic potentials so that two- 
body data still have some significance even for the effective 
ones. Y/e mentroii in this category the worlc done hy Davies , 
Krieger and iiaronger, Iilestor and co-v/orlcers and Donelly, etc. 

The potential of Davies et. al. (Das 66) is a simple 
velocity-dependent potential of the form given hy Green (see 
equation 2.16) . Their potential acts only in relative S- 
state. The parameters are determined from a fit to nuclear 
matter data hut such a potential is in strong disagreement 
with the tv/o-hody data. In an attempt to improve the agree- 
ment with the two-hody data Krieger and co-workers (Krr 66) 
have introduced a density-dependent term in Davies potential. 
The strength and range parameters of their potential are 
determined from a fit to the two-hody data? while the density- 
dependent term talces care of nuclear matter saturation prop- 
erties. The potential of Nestor et. al. (Her 68) is also of 
the form given hy Green hut contains additional spin-orhit 
and tensor terms. These potentials have been used in the 
calculations of finite nuclei (Das 66: Krr 66: Ner 68) in 
the Hartree-Kock approximation and found to give reasonable 
results. 

Brink and Boeker (Brk 67) have investigated various 
effective interactions of different forms, e.g. a sum of two 
Gaussians, a sui:i of a Gaussian and a delta function, and a 
velocity-dependent interaction. The parameters of each potent- 
ial are adjusted so as to give the experimental binding 



energies of end nuclear natter at about the experimental 

densities. With the resulting forces the hinding energies 
and densities of intermediate nuclei are calculated. But 
no force could reproduce all the experimental data at the 
s,ame time. 

Potential Matrix Elements Obtained Directly 

from Phase Shifts 

In most of the nuclear structure calculations the 
primary input needs matrix elements of nuclear force between 
states of two luxcleons in nuclei. The usual procedure 
described above is to construct a potential model from 
nucleon-nucleon scattering data and then use this potential 
to calculate matrix elements needed for nuclear structure 
studies. In an alternative approach relative matrix elements 
of the two-body potential are obtained directly from experim- 
ental phase shifts data without attempting to establish a 
realistic potential form as an intermediate step.. This method 
has been used by Kallio (Kao 65 ) , Koltun (Kon 67 )‘, Blliott 
and co-workers (Elt 67 , 68a). Elliott et. al. have obtained 
the relative matrix elements of the two-body potential in an, 
oscillator basis by using an auxiliary potential method. 

Their matrix elements are commonly knovm as Sussex matrix 
elements . We shall briefly describe below their prescription 
used to obtain these matrix elements. 

The potential V is regarded as a most general non- 
singular potential consistent with the general conservation 



laws of spin, parity and isospin. It may be written as 


V = 


S 


V 


stl'j 


r, |S !/j>< S'U 


(2.31) 


wbere SW denotes the particular spin and angular momentum 
state of the pair and n the oscillator radial quantum number 
describing the relative radial v;ave function R^p(r) of the 
pair. The .Sij labels are referred to as a particular 
channel. Thus, the set of the matrix elements 
<n' Sl!jfVlnStJ> are 




S t Uj 


<n^S L'j!vinS.eJ> = (r 


■^r 


) R^^r dr . 

(2.32) 


In a particular channel V (the label S^J can be dropped for 

brevity if one is confined to a chosen channel) written 

as the sum of tvro terms V = where is the auxiliary 

potential and is small in the sense that it can be treated 

r ■ 9 9 

in Born approximation at some energy E = — , where k is the 

relative momentum and m is the nucleon mass. If 6 is the 
observed phase shift and 6^ that calculated with the auxil- 
iary potential at the same energy B, then it is shown by 
Morse and Reshbach (Moe 53) that 



where is the radial wave ftinction for scattering hy the 

auxiliary potential with the usual asymptotic normalization 

sin (kr - ^1% + 6^) . (2.54) 

The auxiliary potential is taken to he a cut-off oscillator 
2 

^0 “ m" 2:.<a 

= 0 in r >a (2.35) 


with depth D, r.ange a and shape or range parameter h. The 
radial wave function ^^^(2:) is then a solution of the following 
Schrbdinger equahcion 


(iL + 2 IL, + D + 1^2 _ (r) = 0 (2.36) 

dr^ r dr r'^ 4h^ 


for r< a and is the free nucleon scattering wave function 
for r^a. For energy 

™ rn V l^ h^k^ 

S = (2n+'L- j )~ — 2 - - — - = — - ? i2.37J 

mb m m 

the radial part u^,.(n) is given by 

= k(cos 6Qgj^(kr) - sin d^^Ckr) r , (2.38) 

j^(kr) and ■n^(kr) are the spherical j>essel and leumann functi 
defined on page 621 of reference (hoe 53) s is the radial ] 
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of the hound state wave function of an infinite oscillator 
potential.mj^( r) has the required asymptotic normalization 
given in equation (2.34). Constants B and 6^ are determined 
hy the matching conditions 

® ^ lt(cos 6^ o^(ha) - sin 6^ ri^(ka)) , (2.39) 


E^^(a) (cos 6^ 3^{ka) - sin 6^ ii[(ka)) 

“ Tcos““6““l‘‘g^"(ka) - sin 6^ Pg^(ka)) ’ (2.40) 

where the prime denotes differentiation with respect to r. 
Oomhining equations (2.33) and (2.38) v;e have 



r dr = - 


2 

hfk 

2 

mB"^ 




+ ^ Uj^(r)^r^dr. 

(2.41) 

The second term on the right-hand- side ,of equation (2.41) is 
referred to as long-range correction. It arises due to the 
fact that Uj^(r) is not an o.scillator function beyond r>a. 
However, this is very small if a is chosen to be near the 
range of nuclear forces and can be neglected. The matrix 
elements of the total potential 7 are obtained by adding the 
matrix elements of 7^ in the range O^r^^a in equation (2.41). 
In the range a;^_r<oc.>, 7 is same as 7^. Thus, 
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0 -' 








mB 


2 tg(5-o„) 


or more compactly 


V = <BnjJ^o!Sn,>-^'‘E(S-^) • (2.42) 

■’ TTl n 

Tile matrix elements of V are insensitive to the choice of 
auxiliary potential as long as 6 is close to 6^. Equation 
(2,42) gives matrix elements diagonal in radial quantum 
number. The off-diagonal matrix elements are obtained by 
differentiating the diagonal matrix elements vmth respect to 
b. The theory described above has been extended to cover the 
mixing of channels so that the matrix elements in the coupled 
channel can also be calculated. The details can be found in 
reference (Elt 63a) . 

Thus, the auxiliary potential method described above 
consists in assuming that the difference between a rou^ly 
chosen, smooth auxiliary potential, and the real potential may 
be treated in Born approximation. In the earlier work of 
Elliott et. al. (Elt 67) and in Kallio method it is assumed 
that the entire potential V is small enou^ to be treated in 
Born approximation. Thus, the auxiliary potential =0 
and the radial wave function will be spherical Bessel 



fimc^ions for all r. 


(kr) (2.43) 

Kallio uses an approximate relation between spherical Bessel 
functions and harmonic oscillator radial functions given^by 

= B'k , (2.44) 

;. 2 , 2 • , 

E = . (2.45) 


The above approximation is valid at short distances and 
improves with increasing n. Thus^ Kallio matrix elements 
are given by 


<EntlTK,> = 



k 

B 


( 2 


tg6 


( 2 . 46 ) 


Kallio technique breaks down at 6 = %/2 since the matrix 
elements become infinite. In fact Kallio method can be 
regarded as the limiting case of the auxiliary potential 
method described above. 


The Sussex matrix elements are obtained usin.g 
experimental phase shifts of Breit et, al, (Brt 62) and are 
tabulated in reference (Elt 68a) . The range of the oscillator 
length parameter is l,4^b^,2,6B, At each value of b, the 
complete matrix, i.e. for all (n, r/ ) would define the 
potential uniquely. A connection between the n value of a 
diagonal matrix eloment and the energy of the phase shift, 
which governs its value, is defined in equation (2.57). 
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Oonsequently , an upper limil; of 300 MeT, up to wliich experim- 
enteil phase shifts are reliably knovm, iiaposes an upper limit 
on n. Moreover, since the accuracy of the successive dif- 
ferentiation on a curve determined numerically over a finite 
range goes on decreasing, there is a maximum limit of the 
size of the difference n^-n, for which reliable matrix elem- 
ents can be deduced. However, the matrix elements which can 
be obtained reliably are the ones which are of importance in 
most of the nuclear structure calculations. 

In obtaining the Sussex matrix elements it was 
assumed that ■ V is non-singular. However, the matrix 
elements calculated by the auxiliary potential method give 
a first-order approximation to the reaction matrix elements 
even when the true potential is singular. Ihis has been 
shown explicitly by Mavromatis et.'al. (Mas 69). These 
authors start with a reaction matrix (t-matrix) associated 
with a two-body interaction Y and defined as 

t = Y + Y § t . (2.47) 

The energy denominator e is , 

e = Ho - E , (2.48) 

Hq is the unperturbed Hamiltonian of relative motion 

IIq = T + U (2.49) 

with eigenstate 9 and B is energy of the perturbed state 



4 ^ vdiich cp approx 

earlier. If the 

imates. Q is the 

normalization is 

Pauli operator as defined 

chosen such that 

<Ti-> 

= <9j9> » 

( 2 . 50 ) 

then t 9 

= YyjJ 

(2.51) 

and <(9 ft jcp^ 


( 2 . 52 ) 


The matrix elements IT’S ^remain finite even when V has 

a hard core. Of interest for nuclear structure calculations 
are those t-matrix elements for which is the harmonic 
oscillator Hamiltonian, 

= T + i mO^r^ (2.55) 

and the basis function 9 is the complete set of harmonic 
oscillator wave functions. Such a t-operator is defined as 
t^. In the auxiliary potential method of obtaining reaction 
matrix elements, one divides the .potential V into two parts, 

A 

the low-energy part V for which the exclusion principle is 
significant, and the hi^i-energy part (Y-V) , For hi^-energy 
part Q is approximated by Q^, the appropriate operator for the 
two particle system. This is justified since intermediate 
states of high-energy are largely unfilled in the nucleus and 
so exact treatment of Q is not very important. One defines’^,, 
the reaction matrix associated with the hi^-energy part as 
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with = T + V' - E . 


(2.55) 


Keeping terms to second-order in tj^ and Y the perturbative 
expansion (Broii 54) of t^ is 


n 




+ 


± 


A 

V + 



with e 

n 


1 . . 2^.2 
+ niu-J r 


E 


( 2 . 56 ) 

(2.57) 


In Elliott method, Y = Y^ is tahen to he a cut-off harmonic 
oscillator defined in equation (2,55). Ihe basis functions 
cp associated with t^ have the radial part u^(r) defined in 
equation (2,38), If the ' radial part of the wave function*^ 
is denoted 'as W|^(r), the matrix elements 
<Uj^(r) it^(YQ) [Uj.(r)^ = <u^(r) j (Y-Y^) {wj^(r)^ 
may be obtained from the solution of follov/ing Schrodinger 
equations for Ui_^(r) and W]^(r): 




2 2 ^ 

Ju-- + 'll- . ?6 .(- fa ll + Y _ E] ( ruj^( r) ) 
dr El r 




d^ ^ Ktfl) 
dr^ ^ m r2 


+■ Y 


E] ( rwjj( r) ) = 


0 


( 2 . 58 ) . 
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One obtains 

“ j <ir ru]j(r) iV-Y^) M]j(r) 


2 

= - |- k tg(5-5„) (2.59) 

Tlius, using equations (2.38), (2.56), (2,59)and neglecting • 
any long-range corrections we have, to first-order in t/, and 

A 


! "^^nl ^n?^ ^ (6-6q) + ^oI ®n{I> 

+ higher-order terms. (2.60) 

With a suitable choice of V^, it is generally possible 
to make 6^=6 for any particular value of k in which case the 
first term in Gq.uation (2.60) vanishes and I is 
given in first order by ’'^q } even ishen the true 

potential is singular. This is precisely the prescription 
used in the auxiliary potential method -sdiere, however, it is 
assumed that V is well behaved and the two first terms on the 
right-hand-side of equation (2.60) together theh give the 
potential matrix elements ^ * Thus, in any case 

the Sussex matrix elements can be regarded as a first-order 
approximation to nuclear reaction matrix elements , 
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There is a, direoi analogy heiv'/eeiL the Moszhowski- 

Scott separation nethod and the auxiliary potential method 

of calculating reaction matrix elements. In the separation 

method Y is divided into a short-range part V® and the long- 

range part such that the matrix elements of the reaction 

matrix t associated v/ith V® vanish and in first-order t is 
s n 

replaced hy Y^. In the auxiliary potential method Y® can he 

A 9 A 

associated with Y-Y and Y^ with Y, Matrix elements of t, 

A 

associated with (Y-f') can be made equal to zero in which 
case t^ is replaced by Y. However, the separation in the 
auxiliary potential method is not spatial but is energy- 
separation. Unlike the separation method the auxiliary 
potential method can be applied equally well whether the 
potential Y is attractive or repulsive, 

A. 

In Kallio method Y = 0 so that 

• <2.61) 

This is also the expression to which right-hand- side of 
equation (2.60) reduces in the limit 6^ '^0. In fact Sussex 
matrix elements with 6^^ 0 are substantially in agreement 
with Kallio matrix elements. Koltun phase shift method 
is a calculation of t-matrix for the free two, nucleon system 
in an oscillator basis, but we shall not go into the details 
of his method. The work done by Elliott and co-workers is 
most intensive and complete. The matrix elements have been 


tested in some of the nuclear structure calculations for 
finite nuclei and found to give reasonable results. In the 
present worh these matrix elements have been employed in 
the HF calculations of light nuclei in the range 4^A^0. 

II. 4 Summary 

lo summariae, the situation is as follows, She 
experimental information about the nucleon-nucleon potential 
comes from the two-body scattering and two-nucleon bound 
state data. Most of the potential models nowadays use the 
results of meson theory, mainly to determine the long-range 
part of the potential. The short-range part is left to 
phenomenology and the uncertainty in it is mainly due to the 
necessity of worhiiig below an upper energy limit of 300 MeV. ! 
The Ovdgustable parameters of the model are determined from 
a fit to two-body ■ data. One important characteristic of 
these potentials vmich fit the two-body data, is that they ^ 
do not fit the nuclear matter data, e.g. in the binding 
energy calculation of nuclear matter Yale and Hamada-Johnston 
potentials give (Brr 62: Ray 63) a minimum of -8.3 Me'f per j 
particle at a density corresponding to Rormi momentum ■ 
kj, = 1.191'*'^ and -7.8 MeV per particle at = l.llR respect- 
ively, compared to the experimental values of B.E, per : 

particle = -15.75 leV, kp = 1.5R“^. According to Brueckner | 
theoiy, following features of the nucleon-nucleon potential 
are responsible for the smaller values of energy and density: 



(i) large core-radius 

(ii) strong odd-state repulsion 

(iii) quadratic-spin-orlit term 

(iv) weak even triplet central force. 

However, these are just the same features necessary for a 
good fit to hi^-energy scattering data, lahakin potential 
also gives only -8 MeV per particle, in nuclear 

matter calculations, in first order; These potentials when 
used in any perturbation theory treatment of finite nuclei 
fail to give adeq.uate binding in first-order (Ken 66: Shn 67b) 
It has been found that second-order corrections for these 
potentials are q.uite important and reasonable values of 
binding energies of finite nuclei and infinite nuclear 
matter can be obtained only when the contributions of second- 
order terms are taken into account (Ton 64; Pal 66), 

On the other hand there are some effective potent- 
ials, available for use in .binding energy calculations, which 
are computationally easy to handle. The most important 
criteria in choosing these potentials are that they fit the 
nuclear matter data and the second-order corrections to the 
binding energy of nuclear matter are small. Binding energy 
results obtained v/ith these potentials are, in some cases, 
better than those obtained with realistic potentials 
(Das 66: Krr 66) at least in first— order, but then the fit 
to the two-body data is poor. Thus, the question still 
remains open whether a potential model can be constructed 
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wliich. would fit tlie scattering data as well as the nuclear 
matter and finite nucleus data at one and the same time. 

The fact, that different realistic interactions 
which fit the tvro-hody data q.uantitatively give similar 
results in nuclear structure calculations (Lyh 67) » suggests 
that one can avoid the complicated procedure of constructing 
a potential model and obtaining the parameters from a fit to 
two-hody data. Instead, one can obtain the matrix elements 
of nucleon-nucleon potential in some suitable basis directly 
from experimental phase shift data as done by Elliott et. al. 
and others. However, the results of such calculations would 
be expected to be close to those made writh realistic potent- 
ials (see Chapter IV). 



CHAPTER III 


HARTEIE-POCK THEORY EOR EIHITS MJGLEI 


III.l Hartree and Hartree-Pock Equations 

Present microscopic theories of nuclei are hased on 
two postulates 2 

(i) The degrees of freedom associated with the meson 

fields in, nuclei may be replaced by the potentials 
(including exchange forces) between nucleons. 

(ii) Two-body forces are most important. Thus, the 

Hamiltonian of the nucleus consisting of A nucleons 
is of the form 


H = J 


h +> 


i=l 


i<k 


ik 


(3.1) 


Pi 

where T^^ is the kinetic energy of the i nucleon, 
is the interaction potential between the pair of particles 
i and k, and m^ is the nucleon mass. To imder stand the 
nuclear structure one has to solve the many-nucleon Sohrodinger 
equation 


H'JC^, ^2 = E^(?£, rgj > (3.2), 

where ... ^ define the coordinates of the particles. 

This is a non-separable differential equation in 3A independ- 
ent variables and needs some simplification for its practical 



treatment. Hartree and Hartree-Fock (HF) approximations 
(Sir 60) consist of obtaining approximate solutions of 
equation (3.2) using variational principle. In Hartree 
method, the approxima.te gro^md state solution 'pg of the 
above equation is assumed to be given by the simple product 
of A single particle wave functions (orbitals) 

4^2^^) '***^A^^A^’ first A states being regarded as 

occupied so that, 

'4>jj(ri, ... r^) = vL'2(2:2) * 

(3.3) 

The set of functions varied to minimize the average 

energy (pp, system. The minimum cond- 
ition H'^pT P = 0 leads to the following Hartree 

equation for the single particle wave function 

-/■ 




y^aj 


a’ ^1 ^a’ 


d; (?■') vb (x^'')')d? dr" dr'^dr'' 
T a' “ Y a y 


= e 


.H 

a 


Y a Y 

(3.4)..' 


where the summation y runs over only the occupied single 
particle states. It can be seen that Hartree product function 
does not satisfy the Pauli principle according to which 
the wave function must be antisymmetric under the exchange 
of any two nucleons. In the HF approximation many-body 
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v/ave fimction'^ is taken to be fully anti symmetrized product 
of an 0 r til o normal set of single particle v/ave functions | , 

^±(^2) 




^ 1 

'•fir 


, ^ ' “J 










(3.5) 


'I^Js are varied to minimize the energy B =<^^, of the 

system maintaining the subsidiary conditions, namely, the 
normalization of all and the orthogonality of any two 

l.e. 

r ^ _ 

(3.6) 


tp(?) <l7= . 




Thus, the HB eqm.tions are gotten by demanding 

5 [<q 3 , II 4} >/<(];, = 0 . (3.7) 

The details of the variational calculation are not presented 
here and can be found in any book on many-body theory 
(Len 64; Bron 67) . Ihe result is that single particle wave 
functions satisfy the following Hf equations. 

A 

<a|Tjp^ +4) ^ ’ 

where the anti symmetrized matrix element of V is 





AUlP“> = <“Yp|p6> - <^ciy|v-|6P^ 


and <a|I|p> IitpCXi) Si , 

<“y|V|P 5 > = ((4.^(5^) 7 ( 5 J, 4 (?^') 

J P 

dridr2dr:[d?2'' • (3.8) 

The HP self-consistent single particle potential is defined as 


A 


= ^<aY|Tja,jPY^ 

Y=1 


(5.9) 


Thus , 


<d|h|P> „ <<x|T|P^ + <dplP> =s^Ap ’ 

(3.10) 


where h is the rIP single particle Haniiltonian. It can be 
shown that the* u-noccupied HjP single particle orbitals, which 
correspond to the states ay A, p > A, also satisfy equations 
similar to (3.10). That is, 

<a}h!p> = for a>A, p> A , 

but <^a|h|p^ =0/ if only one of the two states a and p are 
occupied and the other one is unoccupied. Thus-, the EP 
Hamiltonian h has no matrix elements betv/een occupied and 
unoccupied levels. The grotind state energy of the system 
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is tlien given by 

K M 

V=1 Y=1 

A' 

= > [<-,il|Y> + I <y|U|y>] . (%]1) 

Y=1 

Equation (3.10) looks like a one-body Sclirodinger equation 
but the potentie.1 U depends upon the very single particle 
wave functions one ha,s set out to find. To carry the sum- 
mation over Y equation (3*9) we need to know the wave 
functions of all the states which are occupied in(p . HencSf 
the solution of HE single particle equ3,tion is obtained by 
the method of self-consistency (iterations) over the wave 
functions. These equations can be solved by choosing any 
one of the following two representations for wave functions 
q^^s, described in next section. 

111*2 Different Yfays of Doing HE 

HE equations can be solved in two different repres- 
entations. 

III. 2a Coordinate-Space Representation 

In the coordinate-space representation the single 
particle wave functions are just written as functions of 
equations (3.8) and (3.9) lead to the following Schrodinger 
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equation when Y is local 


2 

Pi 


A 




Y=---' 


A 




Y=1 -' 


( 3 . 12 ) 


The first term is a local one-hody potential (the direct, or 
Hartree potential) , whereas the second yields the exchange 
term, a non-local one-hody potential (exchange potential). 

If the two-body interaction is already non-local, all one- 
body potentials will also be non-local. 


Calculations in this representa.tion have been done 
by Kerman and Locicett (Ken 64), and Yatitherin and Yene^'roni 
(Van 67). However, they have used a purely central two-body 
potential. Inclusion of a spin-orbit term or a non-locality 
in two-body potential makes the calculations enormously 
complicated* This representation has also been used in all 
calculations of Brueckner et, al. (Brr 61) and Kohler (Kor 65). 

Ill, 2b Oscillator Representation 

The HP single particle wave functions are expanded 
in a convenient set of orthonormal functions 


f 


(3.13) 



where a and denote the angular raoraentun and isospin quantum 
numhers of single particle states, i,e. a = (aljmT), 
(p=n|,jmT): a, n a,re radial quantxam nximhers. Oomhining equ- 
ations (3.13) and ( 3 . 10 ), we get 

y[ <1^112 K> + < oP =eX 

or ^ C“- = 0“ , (3.14) 

where => (Vul 


and 


P. 


\ p 2^*'2 


A 


0 ^ c^, 

1^2 ^2 


(5.16) 


^ is the density of single particle l.evels and the summation 
runs over the occupied states only. The total energy Eq is 
then given by, 


E 


0 



+ i ] . 


(3.17) 


ry 

The expansion coefficients G ’s can be found by diagonalizing 

M' 

a matrix whose elements are . Thus, this is also 

known as the matrix method (Bar 63 ; Net 63 ) of solving HP 
equations. The choice of the basis 9 in equation (3.13) 
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is completely aroitrary “but the choice of three-dimensional 
harmonic oscillator wave functions is preferred for two 
reasons; 

(i) The HF single particle functions are expected to he 
close to oscillgitor ones, so that a few terms of the series 
may be needed to get a good representation. In other words, 
the set of oscillator functions used need not have a large 
dimensionality and instead of solving integro-differential 
equations of the type (5-12) one has to diagonalize small 
real symmetric matrices. Thus, the calculations are much 
ea.o.ier. The coefficients and energies e’s are always 
real since they are the eigenvalues and eigenvectors of a 
real symmetric matrix. 

(ii) In this representation the sepa,ration between the 
relative and centre-of-mass coordinates can be done exactly. 
Thus, the calculation of ty/o-body matrix elements of nucleon- 
nucleon potential is easy since the latter generally depends 
only on the relative coordinates of two nucleons. 


Tbis representation has been used in most of the 
earlier HF calculations described in Chapter I as well as in 
the present work. The function o is the normalized solution 
of the three-dimensional isotropic harmonic oscillator; 




m 

eigen function of total angular momentum defined 

in (2.14), 7^ is an isospin function and the radial function 

is 




2fp^(n-l)i> 




p(n+e^-)^ J 


cm 


yi^,+i ^-x/2 mm 


■'n-1 


(X) 


(3.19) 


2 

X = p^r , n = 1, 2, 3 Ii“(X) is a Laguerre Polynomial, 

defined on page 784 of reference (Moe 53) and 
being the oscillator frequency, ' Ihe energy eigenvalue 
associated with this state is e^^.. (2n+t-^)hui. 


In general, the summation over p in equation (3.13) 
is over all the quantum numbers. However, the HP Hamiltonian 
h can have various symmetries, such as spherical, ellipsoidal, 
axialj'etc. v^hich when considered in the initial step, can 
reduce the vast vamiational space in expansion (3.13). If 
we consider the HP field to be spherically symmetric, then 
the single particle HP states have good angular momentum and 
the 0^ (and hence A, ,,/) are diagonal ±n t , j and independent 

jj, ' ^ 2^2 

of m (z-projection of f) . Thus, the expansion (3.13) can be 
rev/ritten as; 

. H 

4^a(r)' =^ae3mT(^) = C!|( tjv) * (3.20) 

where H specifies the maximum number of terms which are 



included in the suLoiiation in equation (3.13)» and is same as 
the maximum nuciher of relative nodes in functions occurring 
in the summation, Ihe mixing of neutron a,nd proton states in 
a given orhit has not been considered so that expansion of 
equation (3.20) is limited to a single value of third comp- 
onent of isospin, t which distinguishes the neutron and proton 
orbits. Also, if there is no coulomb force G^s are indep- 
endent of T . 

Ill . 3 Calcu l atio n of the Matrix Elements 
in 0 sc iiratb~r Basis 

III. 3a One-Body Matrix Elements, } '^ • 

The matrix elements of the kinetic energy operator 
T in the oscillator basis is given by 

6(3i3i) 6(m3_m5_) )( 
n^ [ T ( J n^^ , 

where 

^ (2n^+i^— ■§■) , n^=n^ 

f 

! 

■(^n^l T(i 2 _) j ^ yV ,|n^-n^j =1, n=min (n^n^) 

j o * ^5-21) 

III. 3b Two-Body Matrix Elements, * 

Since most of the two-body interactions depend 
only on the relative coordinates of the two nucleons, it 
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is convenient to 
terms of relativ 
tvro particles.. I 
■body uncoupled st 


erqpress the two-body matrix elements in 
e sjid cent re- of -mass coordinates of the 
he follov/ing procedure Ib used, The two- 

--.■o-e s I ^ ! ^1 3 l^'l'^ 1 » ^2 ^2 ^ 2^2"^ 2 y 


are coupled to a total angular momentum J and total iso spin 
■ 1 state with tiie use of 3-0 symbols; 




> (- 1 ) 


2 “rJ 4*Ll4'i-'.v 


Ij]* [r]» X 


JM'CM 




/ X 32 j\ /» 4- <\ 

y ! \ i 

' - 1 -m^ -m. M 1 It. t. -11. 


‘1 ^2 




\ 


A 


1 0^ » n2^/20 2” 1®^? h > 


( 5 . 22 ) 


vfnere we have used the notation that a number x inside sq.uare 
brackets implies 

[x] ( 2x+l) . 

The two— particle couj)led state in o^O coupling scheme is 
transformed to the >,-3 coupling scheme by means of 9-0 symbols 
in the following wayt 
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X • (3.23) 

« 

Further, Moshinslvy transformation is used to write 


~ ^ nLNl!^|n^^2^n2i2?'X I XX » 

nfclTL 


(3.24) 

where the n^IL ?\jn^|,^n 2 |, 2 )\^ 3-^® Brody and Moshinsky 
transformation hracliets (Bry 60) . The quantum numbers (n|j) 
and (WL) specify the relative motion of the two particles 
and the motion of the centre-of-mass of the pair respectively* 
The total orbita.1 angular momentum the result of the 

coupling of £ and L. A simple recoupling procedure yields 


|nt, NL, i t, S: JM)> - 

= (-1)^+^^^^ ([MLj])^ ^ X nisiX jm) . 

J' - Is J JM 

(3.25) 

The quantity inside the curly-b rackets is the Wigner 6-3 
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symlDol, The final expression for the most general antisym- 
metrized tv/o-hody matrix elements of an interaction V which 
« 

conserves total spin S and iso spin T" is 


^^ 1 ^ 1 ^ 


l^l'^l^l''!^ 


^2^-J2^2^2^2\ ^2^ 2'^ 2^2^ 2)^ 




Mi:’ 

I^nn't IJ m 


X<ntlTL Xln^^ 1^2^ fy ^ ^ ^ ^1 . X 





f 


( 3 . 26 ) 



where 


J S * ' 

T \Y\n' ^ Si J M,^''> 




( 


^nlSJlfS 


f2 


) V(r) gjj^( 


— ^ 

r 

f2 


)dr 


\h> 


(3.27) 


In most, of the situations, the interaction V(r) conserves 
total angular momentum, parity and isospin. The 3-2 t 6-3 and 
9-3 symhols follow the notation of Edmonds (Eds 57) . The 
normalization constant for the matrix elements is 


(1+67 n ^ 1/ i 
^ 1^2 ^ 1^2 




Thus, apart from the geometrical factors occurring in the 
expression of equation (3*26), the two-body matrix elements 
can he obtained in a straightforward manner from the relative 
matrix elements of V in an oscillator basis. 


111.3c Matrix-Elements of HE Single Particle 
P 0 tentia,! , ^ ^j_i ^ | ^ * 


Since the 


summation in equation (3.20) is restricted 
to the radial quantum numbers, only the matrix elements of 

the form<n^73ii'^l'C^, n2?,232^2'^2| ^2^232“2'^2> 
are required in the calculation of • Y/riting 


explicitly, 







2Z. f 

*^2=1 

( occ) 




^i\22) X 




( 3 . 28 ) 


with. 


° ^ “n,(t232) <’S^2h^ ' 


a( occ) 


In the case of a douhly closed shell nucleus all the occupied 
j 2 “ suId shells are completely filled, i.e., each one of them 
contains 2(2o2+l) particles, Thus, the summation over m 2 and 
T 2 equation (3.2$) can he carried out very easily and the 
orthogonality of the 3-D symbols loads to the result 


/ n(t<232'”2^2> 


^ 2 ''' 2 


^^^1^^ SCDqDq) (2D2+1) ^^^1^1"\^1* ^2^2^2^2^ 

( 3 . 30 ) 



b!? 


nn'i ?,' HL 


(-)''+ ^'[J'lLAlL XlLSlEJlLr] /n&UMn-LC^ng?,^ 




1^2 

M 

i 

X <(n' [! hL A \ ^-[1-^2 ^'2 ^ y 

)* 

1 

1 

2 

t 

s V 

^1 

1 


j 3l 

jg 

J i 

J 

V* / ~\ ( 1 'i ^ +S-i^'^'\ ( 1 J s 

A ^i-(-i) V^imat^ 

i L 

i 

J 

! 

J S 

J 

t/ ?" 

J ! 


j 

4^2 

a '., 

1 

] 

) i 

1 ^ 

1 

2 



m 

■ ^1 

V 

^2 

\ 


u t' >' 'i 
Is J jM 

V. 

(3.31) 


with 


pJ^T s 
Wl- V 




S 

nn' t i! 





1 


J^!7~ g 

nn' L I ' 


(3.32) 


Thus, 


= 5(di2_^^) 

<ni|U(!,iDi)|ni> , 


(3.33) 
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where 






Pn^n'^'-kh'i ( 232 + 1 ) 


P ( f •1^2 3 2^2^2^ * 

(3.34) 

For non-closed-sliell nuclei, which contain partially filled 
j-suhshells, following approximation is used for the density 
of the occupied levels. It is assumed that the neutrons 
and protons are equivalent and averaging over the m-suh states 
for a given j ic done. The average density is given hy 




= r 


a( occ) 


r^a. (f . ^ | - x 

2( 23^+1)' °n„^^232^ ^n^^^2^2^ 


(3.35) 


where h (loOo) = number of nucleons in a given suh- 

3/ 41 c 

sIigH ^^2^2^ * 


Equations (3.29) to (3-33) remain the same. 


III. 4 


Iterative Solutions of HE 



With equations (3.21) and (3*33) s equation (3.1-4) 


becomes 




where N is the n-unber of terms included in expansion (5- 20). 
There is one such equation for each . The total ground 

state energy defined in equation (3.17) is 




B, 


71 -I 


n. 




i n-,n:( ^^1^1^ 2(2j-,+l) \ 
* 1 1 


/ 






(3.37) 

The factor 2 ( 23 ^+!) comes from the sum. on m^ and The set 

of equations (3.36) is solved by the following iteration 
prooedure : 

(i) An’ initial set of guessed and also which 

orbits are occupied. Usually, = 6^^^. 

(ii) With this set of coefficients the matrix elements 
of HE Hamiltonian occurring in equation (3.36) are 
calculated. 

(iii) Equation (3.36) is solved by diagonalizing a IXH 
Hamiltonian matrix. A new set of coefficients 

obtained. Total ground state energy 

is computed. 
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These new set of coefficients are recycled throng the step 
(ii) for recalculating the matrix elements and another diag- 
onalization of the matrix. This process is repeated until 
successive diagonalizations produce the same set of coef- 
ficients and reaches a constant value. The final 

set of correct expansion coefficients} 

is the total binding energy of the system and 
single particle levels of the nucleus. The calculations are 
repeated for several values of oscillator range parameter 
b = Y'(fL/m^U)) , vdiich enters in our expression throu^ the 
expansion (3*20). The minimiim in the E^ vs. b curve is 
determined. 

Ooulohb and Centre-of-Mass Correction s 

There are two important corrections to HE solutions 
arising due to the coulomb force between the protons and 
centre— of-mass motion of the. nucleus. The effect of these 
corrections on the total binding energy can be estimated 
externally if one assumes some simple model for the nucleus. 
Thus, for example, if the nucleus is assumed to have a uniform 
spherical charge distribution, its coulomb energy would be 
given by 

2 = 3 ^ Z(Z-l) with R = 1.3A^^^ * (3.38) 

•^coul 5 R 

Similarly, if one assumes that the centre-of-mass of the 
nucleus is in the lowest Is state of a harmonic oscillator 



well, the kinetic energy of its motion is 

( 1 - 59 ) 

However, the effects of these corrections on various 
nuclear properties calculated in the framework of HF theory 
can he determined only if they are accounted for in a self- 
consistent manner, i,e., these corrections are included in 
the Hamiltonian h first and then the equations are solved. 
Ihr this, the following procedure is used. 

III. 5a Goulomh Correction 

The Goulomh force between the nucleons can he 
wristen as 

. ^ 5 . 40 ), 

^coul particles are protons. The Hamiltonian 

.H wiljli "tliB couloml corr6c1jioii is 


H 


JX 


Tj_ + 




i<_k 


r 

i=l 


Ti + 


e: 

i<k 


ew 


(3.41) 


new _ Y + ‘ 

where " vik ik 


(3.42) 
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Tlius, the only change which would oooui In the HP equations 
is that the two-body matrix elements <hit‘2l’^lk|‘*l^‘2> 
will be replaced by those of vJJ". These matrix elements 
can still be given by expression (3.26), but now will 

XlXl V L- 

consist of two terras, 




,J ' C 

nn'’ t (f 


S _ S 


(nucleus) + S 


nn'' It 


(coulomb) ( 3 . 43 ) 


with 


J^'T S 

\n'’jti' (coulomb) = (ntsJM, TM ‘ ® 


T ! r 

QG 

( 

0 . 


r/l' SJM, TM^> 
2 


= 8('n) &(M^)5(|,£,') f- 

(3.44) 


Care should be excercised in carrying the summation over ^2 
in equation (3-30) since coulomb matrix elements depend upon 
the isospin quantum numbers, Hovrever, the sum over Tg can 
still be carried over explicitly if v/e malie the time-saving, 
but non-essential approximation that in calculating , 

P /(todo) independent of Vp- Summation over mp is dona as 
before, Thus, equation (3-31) will now be replaced 

by 

'7' ^ =r (nucleus) + (coulomb) (3*45) 

rm'Ci: ' im IL 
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with. 


gj'r s 

^nn'/. f/ 


(coulomb) =2 





2 





j'T'S 
nn't t' 


S(Tni) S'! 


./t s 


(3.46) 


Subsequent changes will occur in the expressions that follow 
equation (3,30), 


III. 5b Gentre-of-Mass Motion 


Considering a nucleus with A particles, we use 3A 
coordinates for A orbitals-. However, the centre-of-mass 
coordinates should be subtracted, so that we really have 3A-3 
independent coordinates. Thus, the Hamiltonian Miich desciibes 
the intrinsic state of the system is obtained by removing the 
centre-of-mass energy in equation (3*1) • Thus, the intrinsic 
Hamiltonian is 



(3.47) 


with the total momentum ^ • (3.48) 

i 

The operator of the centre— of— mass kinetic energy can be 
written as a sum of one-body and two-body operators. 


' p2 

2m. A “ fL. 2m.jr 

^ ik ^ 
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The square of the relative momentum is 

^Pik^^ = I'Pi 

= - 2pT.^) 

••M*^k = K'-Pk) - pL 

^^Pi+Pfc^ - ^ * (3.49) 


Equation (3.47) novf becomes 


Hx = 


1 ) 

2 i‘VS: 


''ik 


-I 


1 . 


gpL 

-to. 


(3.50) 


which is a pure two-body operator. The secondpterm can be 

p 2p-| p 

regarded as an additional potential V^p * = and ©an 

k. 

be treated as the nucleon-nucleon two-body interaction 

g 

This adds an additional term to in equations (3.26) 

and (3.43) j which is now given by 

(o.m.) = 6(U,7f <n|l(t) In') , (3.51) 

where <(n|T(t)tn') is defined in equation (3,21). The solution 
proceeds as before except that now in (3.36) and (3.37) one- 
body term <(nj T(£^) |n^ is missing. 



III. 6 Hartree-fock Density Distribution 
and R.M.S. Hadii 

III .'6 a Density 

Ihe density of the particles inside the nucleus 
is 


? — ^ ^ 






•*• ^^a-l'^^a+l °* 


(3.52) 


Using exp 


cpression (3*5) fo^ Q3 integrct 


ting over 


(A-1) 


)i:diiia'teS| one olDtains 






a = 1 


= > P,,' • 


(3.53) 


Writing and cp^Js expl- 

(3.18), and summing over m 


ss explicitly from equations (3.16) and 


and T we odisain 


e4^(3^) 

u)=fei: 

Wtj 






>.54) 


This density i£ 


normalized so that |r)dr= A. The proton 
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and neutron charge densities, f„(r) and P(r), have the same 
form except that in equation ( 3 - 53 ) the summation is only 
over the proton or neutron orbits respectively, and the factor 

ecuiation ( 3 . 54 ) is replaced by ) • 

Finally, the three radial densities are related as 

f„(r) = p^(r) + f^(r) . (3-55) 


111 . 6 b R.M.S. Radius 

Ihe mesn square radius for the matter distribution 
is given by 


<h> = 


II 

2 

r dr 

(3.56) 




\ * 



= (p 1 

A /\ 0 

(r) r"^ dr • 

(3.57) 


Using equation (3.54) this becomes 


y A (2o+l) Pnn''^^^) 1 ^ 


nn^'i'^ 0 

(3.58) 

The integration can be carried over exactly and finally 
(r^)^ ^4b_y— ( 23 + 1 ) [^^(t^) ( 2 n+f^|-) 




“ n+ 1 ^^ 3 ) fn(n+&4%) . 


(3.59) 
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2 

Tile r.m.s, radius is '{ <^r ^ . The mean square radius for 
charge distribution is given by expression similar to equation 
( 3 . 59 ) except that A is replaced by Z and sum is only over the 
proton orbits. 

111*7 Eg Potent i al 


The oscillator matrix elements of the HE single 
particle potential U are obtained from equation (3*34) after 
several iterations have determined a self-consistent density 
matrix is possible to express U in coordinate 

space by the equation, 


.11 


TD'(r|rO = 


j 1 

^1^1“^ 





(r) 


< 






( occ) 




Using explicit forms for 9 ’s we have 


U(rl?0 


‘'1 “1 


biDimiXi 


^1^1 




% 


2 


(3t6l) 


where 




= 1: 




-i-i — — 5? — 

(3.62) 
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This is a non-local and state-dependent potential. A possible 
local eq_‘u.ivalent to tnis potential can "be obtained in the 
static approximation which consists of computing the velocity- 
dependent potential U(rj^ corresponding to non-local U at 
^ = 0 . Thus , 

U(r'j^ = dr" 



(3.63) 


after doing the angular integrals. The integral on r can be 
done by using the Laplace transform on page 30 of Magnus and 
Oberhettinger (Mas 49) . The result is 


r 

0 ^ 




rdr = 


(n-1) ! 



(3.64) 


The static potential becomes. 


u(?!?= 0) = (-) 


n 


1 






< n^. |U (Oi) I n^) b5/2(|) V4x 


n-, n 





(3.65) 
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III .8 Choice of lltiolei for Spherical HP Oaloulations 

The HF‘ Hamiltonian h does not necessarily have the 
same symmetry property as the nuclear Hamiltonian H. The 
possible symmetries of the HF Hamiltonian are governed hy the 
following theorem (Hia 68). 

Theorem: let R be an operator that commutes with H, 

-1 

R HR = H . (3-6S) 

Then, if R leaves the set of occupied orbitals invariant, it 
commutes with HR Hamiltonian h, 

R”^hR = h . (3.67) 

The operator R v;ill leave the set of occupied orbitals invari- 
ant if the states obtained by operating R on any occupied 
orbitals a may be expressed as a linear combination of occupied 
orbitals only. That is, 

A 

0 = 1 


“4 hi 

a = 1 

If the nuclear Hamiltoniaj' H has spherical symmetry, it com- 
mutes with the three operators (Tx* fotal 

angular momentum operator We wish that HR Hamiltonian 
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should also have spherical symmetry in vmich case it should 
commute with and These operators will leave 

the set ia = aijoiTl. of occupied orbitals invariant if for 

j 

a particular 2 all m-suhstates are either filled or empty. 

Hence, spherical sjnmietry can exist only in nuclei that iiave 
right number of particles to form closed shells. Such nuclei 
in the range A ^ 40 are %e, ^^0 and “^^^Ga and, thus, one 

would expect that results of spherical HI? calculations would 
describe the ground state of these nuclei correctly. If the 
coulomb interaction is neglected and nuclear forces are assumed 
to be charge independent then the HE field is an iso spin scalar, 
and neutron and proton orbits are degenerate. 


111*9 Single Oscillator Approximatio n 
Limit to "HE Solutions 


In the limit v/hen there is only one term in the 
expansion (3.20), the HE orbitals are represented by pure 
oscillator functions. In this case P .. .. i will become 


f 


^ 2^2 


V 

z_ 


a( occ) 


,a 


G 


/ 

2 



a( occ) 




2 


a( occ) 


6 6 / • 

ocl-^2 ^2^2 


( 3 . 69 ) 


Thus, the density P, , / of the levels would be either ene bf' 

^ 2^2 '■ ^ ^ 
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zero depending upon wlietlier ilie state pg is occupied or unoc- 
cupied in tile ground state configuration. The single particle 
potential matrin elements are 





^ 2 ( 000 ) 


The total energy HI- is 


P2M-2 


(3.70) 


E, 


f *^lf^j^l> 2 <^“'''l|^|^l > ^ 




1 

2 


p^(occ) 


o 




^il ® 1^1)^. 2 <! *-''if^ 1^1)^ ^ 


(3.71) 


where the potential energy of the nucleus is, 

^ ^ ^ ^1^2 j^AN |^1^2 ( 3 . 72 ) 

Eor If = Z nuclei with IS closed shells equation (3-72) reduces 
to a simple expression. 





= 2^ nSiL <ii^SJjV|ntSJ> , (3.73) 


i — — 

St: J 


where the v/eigiit-f actor , lo given hy 


^iSTJ _ /l-C- 


0 




■) [V] [J] X 


\ . 

[>■1 

V. 



( occ) 

The potential energy can also he v/ritten as 

P.E. = 

p(occ) \x ^ 

( occ) 


(3.74) 


(3.75) 


where W(p, ) is the eontrihution to potential energy due to 
interaction of nucleons in the state p, and can he obtained 
with the help of equation (3.72). for douhly closed shell 
nuclei, the single particle potential U, defined in equation 
(3.70), can he written down esjplicitly as, 

{nq 2 1 i ^ I ^ 2 “ 2'’' 2 ) 

[£.1 LS <;n3_i^ji, n2£,232- M^jY 


S~ Y~ lJJ 

^~~S~ • ■^;=p. 1- ^ X-l 

no to 3 o Jv 


2 ^2 '^2 
(occ) 


[j^ni^iJi, 

+ (-l) 




| n 2 t 232 ’ ^ 1^1 ^ 1 * 

(3.76) 



Equation (3.76) defines th.e single particle potential for 
neutrons. For protons, additional coulonl term should he 
added and the result is that Y in equation (3.76) is replaced 
"by defined in equation (3.42). fhe kinetic energy T of 

a single nucleon in state 1 is given hy 

f = (2n. + - -1) . (3.77) 

h 


III, 10 Validity of the ijE Approxima.tion 

The deterninantal v/ave function constructed from 

the final self-consistent HE single particle orbitals 9^s 
cannot be the true many- body y/ave function and hence it is 
relevant to ask the question hov/ good is Yfith the pos- 
sible improvements to one can y;rite an improved as 


vjj - -^Q + ■^l^^HE’ 


,2p-2h 

HE 


(3.78) 


where 





2p-2h 

HE 


~^^ph ^-p^h ^HE * 


= 2^^p^P2:hih2 "Pi"p2%"h2 "fHE 
^ 1^2 
■h^ig 


( 3 . 79 ) 

(3.80) 


p and h refer to the HE unoccupied and occupied states 

respectively. The various amplitudes, i.e., , A-,1 , . 

O’ 1 ph ’ 



82 


2 P3_P2- 

large compared 


Making use of tlie 


are in principle calculable. If Aq is very 
the other amplitudes, ^ 
fact that the HP Hamiltonian h does not 


close 


Yr 


have any matrir elements between the occupied and unoccupied 
states, it can be shovm that the amplitudes for the Ip-lh 
excitations are identically zero and, thus, the corrections 
to can come only from 2p - 2h and more complicated 

excitations. She admixture of 2p--2h excitations to is 

measured by the ratio of some lypical off-diagonal matrix 
elements of the Hamiltonian H, between the wave functions 
corresponding to these excitations and bo bbe smallest 

amount of energy needed for two particle excitations. A small 


ratio results in small admixture. Same is true for higher 


excitations. Shis, then means, that the validity of the HP 


approximation is measured by the difference in the HP single 
particle energies of the unoccupied and occupied states. She 
lowest of these differences is the one between the uppermost 


occupied level and the lower most unoccupied level. If this 
is significantly large compared to the differences in Idle 
energy of the unoccupied states themselves or of the occupied 
states, we can speak of the existence of a gap of energy. 

She larger this gap is, the nearer is to Sometimes, 

a relaxation of the constraint imposed on single particle 


orbitals results in a considerable improvement in the 


approximation and is reflected in increase in the energy gap 
between occupied and unoccupied orbitals. 
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HP appr-oximation is a variation -.procedure as described 
earlier in Section III.l and provides an upper bound to the 
energy of tiie system. The possible corrections to arise 
from 2p-2h and higher excitations, and can be calculated by means 
of perturbation theory in Hartree-Eociv basis. However, HP sol- 
utions have to be imderstood in a different way vdien one is 
dealing with the effective Hamiltonian of Shahin, Waghmare and 
co-workors. Those authors obtain an effective Hamiltonian by 
using the theory of cluster expansion and higher-order terms 
in this expansion are dropped. HP method to such an effective 
Hamiltonian has not been applied as a variation procedure since 
we are not dealing v;ith complete Hamiltonian. We, then, have 
used the HP method as an attempt to find orbitals in this 
harmonic oscillator basis that yield self-consistency so that 
Ip-lh corrections to the energy, calculated with this effective 
interaction are minimized. 



GHAPTEE IV 


APPLICAEIOIIS OP EAEIEEB-POGK IHBOSY AMD DISGUSSION 

IV. 1 ' Details of ijlie Galoulation 

Spherical Hartree-PocA (HP); calculations have "been 

A R 

carried out in oonsideraole detail for the nuclei ^He, Be, 

12 16 

G and 0 wi'ffi effective Tale interaction. Less complete 

1 0 18 

calculations have .also "been done for “ Bg and 0. In view 
of the fact that Besides the properties such as Binding energy, 
root-mean-square ra,'dius, spin-orbit splittings, etc. , we are 
also interested in the distribution of protons inside a nucleus 
and that the notion of the centre-of-mass of the pair of part- 
icles contributes considerably to the binding energy, we add 
the coulomb and centre-of-mass correction terms to' V vdiile 
calculating the HP potential and wave functions (see Section 
III.5 for details), Ihe effect of these tvro terms modifies 
the wave functions through the density matrix. These effects 

have not been included in a self-consistent way by Shakin 

16 

et. al. (Shn 67b) while calculating the properties of 0 
and ^^Ga with effective Tale interaction. However, while 
they have included the second-order term in pseudopotential 
(VP) by adding a term VP VP to defined in equation 

(2.27), we have not included it in our calculation. HP 
calculations have also been done for the nuclei '^He , ^^0, ^^0 
and "^^Ga with Sussex interaction. 


In these calculations also. 


coulomb and oeiatre-of-mass matrix elements 


are included in 


a self-consistent manner in all the cases, except in 
for which the effect of these terms is estimated externally. 
Iwo-body matrix elements are calculated using formula ( 3 - 2 . 6 ). 
Since the nuclesn force has a short-range it is necessa2:y 
to use only a few relative ^ values in the calculation of 
two-body matrix elements of the nucleon-nucleon interaction. 

A maximum of relative 1^=2 states has been used in our oalcul- 
ation. Higher E states are not expected to contribute signif- 
icantly for light nuclei studied in the present worli (Tan 68). 
The configuration space used in our EH calculations of A ^ 16 
nuclei consists of s^/2» P3/2 /2 states of 

n = 1, 2 and 3 major shells for each of these states. In the 
calculation for "^^^Ca and also for some of the lighter nuclei, 
a larger configtiration space is employed which includes ^^^2 
and oscillator states of n = 1, 2 and 3 major shells in 

addition to states described -above . Thus, in all the HH 
calculations expansion (3.20) contains only three terms. 
Earlier HP calculations (has 66: Sve 65s Shn 67b; Krr 66) 
provide some justification for such a truncation. However, 

this truncation r.iahes our results depend on b( = f )» 

A 

the oscillator range parameter and in the spirit of a proper 
variational calcul-ation the criterion for the best choice of 
,b is the one which minimizes binding energy. After a suf- 
ficient number of iterations have been performed, the computed 
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q.uantities converge to constant values ,and the numher of iter- 
ations needed for such a convergence depend upon the initial 
G-coefficients. In all the calculations r/e start ?/ith initial 

iterations have been found to give a 
satisfactoi^y convergence of computed cuantities which finally 
describe the properties of the nucleus. In order to study the 
nature of the HF wave functions and also the single particle 
potential, we have concentrated mainly on the results obtained 
for 0. Ilie reason for this being, although there is quant- 
itative disagreement between the HF v;ave functions of various 
nuclei studied here, the general nature of these wave functions 
remairs the same for all the nuclei. The discrepancies wher- 
ever encountered have been pointed out. Single oscillator 
configuration approximation to the HF solutions has also been 
discussed. We present the results of our calculations in the 
following sections. 

I'^* 2 Effective Yale Interaction Results 

IV. 2a Binding Energy", Single Particle Energies 
and Root-LIeaja-Square Radius 

The ground state total energ;/- of the nucleus is 

calculated as a function of oscillator range parameter b. The 

results of the binding energy as a function of b are shown in 

A R 19 1 

Figure IV. 1 for the nuclei He, Be, C and 0, Ideally, 
the HF resxfLts should be independent of b or basis chosen and 
Eq vs. b curve should be a straight line. This would be true 
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only when the RF wave functions are expanded in a complete set, 
Ihe larger the number of terms included in expansion (3.20), 
the ¥/eaker would he the dependence on b. Hov/ever, from the 
figure we notice that the curves are reasonably flat near the 
values of b which give maximum binding. She two-bodj’' matrix 
elements needed for these calculations are evaluated using 
tables of relative matrix elements for the effective lale 
potential given in reference (Shn G7b) . 

In Table lY.l, we list calculated HF properties of 
various nuclei for the values of b which give maximum binding 
(b=1.76F for "^He and b=2.09F for ^Be, and ^^0). These 

properties are the binding energy per particle (B.E./a) , 
single particle energies, matter and charge root-mean-square 
(r.m.s.) radius, and spin-orbit splittings (i-s splittings). 

The corresponding experimental values -wherever available 
are also given for the sake of comparison. Experimentally 
single particle neutron levels have been determined from the 
stripping process, as for example, the (d, p) reaction. 

Analysis of the experimental data on cross-sections for 
excitations of the various states of the final nucleus and 
angular distributions corresponding to these states gives 
information about the shell model single particle neutron 
states. Similarly, neutron hole states are. located with 
(p, d) or (d, t) pick-up reactions while proton hole states 
are determined v/ith the help of (p, 2p) and (e, ep) knock-out 



Calciolated HF properties with effective Yale interaction for various nuclei 
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Riou (Riu 65) ■ d. Hofstadter (Hor 67) 

Cohen (Con 63a, 63h) e. Ulrich et. al. (Ulh 58) 

Mattauoh et, al (Mah 65) f. Vaahakidze ct, al, (Vao 68) 
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reactions, Prcti tlie Ta'ble IV. 1 we see that "binding energy 
results are in poor agreement with the experimental values. 

‘The reason for the discrepancy in the binding energies ster^ 
from the fact that we have not taicen into account the second- 
order term in YP. Calculations of Shahin et. al. (Shn 67b) 
show clearly that rhis term is isfHportant as far as the bind- 
ing energies of nriclei are concerned. In fact this term 
alone gives a binding of about 5-4 HeY per particle. In 
view of this, the theoretical results are in satisfactory 

agreement for and ^^0 but still somevdiat poor for ^Be and 
12 

0. This is, hov/'ever, not surprising. Prom the experimental 
information on transition rates and quadrupole moments these 
nuclei are knovn to be deformed. This conclusion is also 
reached from the ca,lculations of Bassicliis et. al. (Bas 67). 

In fact, these authors have made deformed HP calculations 
(mixing different orbital angular momentum states) and obtain 
binding energies for °Be and G which are lower than spher- 
ical HP estimates. Their results are in satisfactory agre- 
ement with the experiment. Pal et, al. (Pal 67) have also 
recently performed similar calculations for the 8^ 1:^40 
nuclei and obtain satisfactory agreement for the deformed 
nuclei in this region. Eo?7ever, these authors have not 
included coulomb md centre-of-mass corrections in a self- 
consistent manner. On the other hand, if we include more 
orbitals in the expansion of basis states our results would 
be somewhat improved. 
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lPiy '2 splitting 
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"better if second-order terns in VP are included (Shn 671)). 

Pile HP single particle levels are in satisfactory agreement 

with the experincntal values. There is e. large gap betv/een 

the occupied and the unoccupied levels in the case of "^He and 
I'G 

0 and is a characteristic of HP approximation. The calcul- 
ated charge and matter r-m.s. radii a,re in fair agreement with 
the experimental data. When the second-order term in VP is 
taken into account the particles interact more strongly and 
it is expected that this v;ould result in sn increase in the 
"binding of HP single particle levels and decrease in the 
values of charge and matter r.m.s, radii. 


In order to estimate the effects of couloiii"b (coul) 
and centre-of-mass (c.m.) corrections on various nuclear prop- 
erties we have repeated the HP calculations for the following 
cases. 

Case A: Without coulom"b and centre-of-mass corrections* 

Case B: With centre-of-mass and without coulom"b corrections* 
Case Gi With couloa"b and vrithout centre-of-mass corrections* 
Case J)t With both the corrections included* 

Table IV. 2 shows the total ground state energy obtained 
in different cases for b=2.09P. The kinetic energy of the 
centre-of-mass motion is determined from the difference 
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Det5.7een obtained in cases D and C. Sinilarly, the coulomb 
energy of the nucleiis is obtained from the difference of total 
energies in cases Jj and B. From the table, we sec that 
centre-of-mass correction is very important for the nuclei 
He, Be while coulomb correction is not so important. For 
coulomb and centro-of-mass effects tend to cancel each 
other. These corrections to the binding energies have also 
been estimated enternally by assuming that the charge distr- 
ibution of the nucleus is that of a uniformly charged sphere 
and centre-of-mass moves in Is state of a. harmonic osc-illator 
well (see expressions 3-38 and 3*39) and compared with the 
values obtained self-consistently . It is seen that in all 
the cases the spherical model overestimates the coulomb energy. 

The effect of these corrections on the HF single 

particle energies is shown in Table I¥.3 for b=2.095'. The 

energy of the occupied state for ‘^He becomes almost 

double ?;hen the centre-of-mass correction is included. This 

correction is ciso very important to reproduce the correct 

12 

ordering of Ip^^o ^^1/2 particle levels in 0. 

In Table IV. 4, we give the matter r.n.s. (r^) and charge 
r.m.s. (r ) radii obtained in different cases for b=2.091. 

The coulomb correction increases r^ in all the nuclei, whereas 
the centre-of-mass correction reduces r^ in the case of ^He 
and by approximately 20 percent and 8' percent respect- 

ively. Besides, charge r.m.s. radius is larger than matter 



Effect of coulomb and oentre-of-mass corrections on HE single particle energies 
of various nuclei obtained with effective Yale interaction, b = 2.09E. 

iU.1 energies are in MeV. 
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TalJl e IY.4 

Effect of coulonlD and centre-of-iaass cc rrections on 
charge (ro) and laass (r ) r.n.s. radius of various 
nuclei obtained v/ith effective Yale 
interaction, b = 2.0SE. 


Nucleus 

Case A 

Tq (E) 

Case B 

I'o 

Case 

rjl?) 

0 

^0 U) 

Case 
ro (E) 

D , , 
(?) 

^Heg 

2.2229 

1.9013 

2.2382 

2.2438 

1.9066 

1.9096 

"Be, 

5.0335 

3.3531 

3.0873 

,3.1113 

3.4283 

3.4735 

12p 

^"6 

2.7403 

2,7526 

2.79 26 

2,8144 

•2 . 8158 

2.8571 

IGq 

^8 

2.5420 

2,515v 

2,5734 

2.5876 

2.5486 

2.5699 


centre-of-niass corrections on HI’ calcnlations has also been 
recently studied by Gunye (Gue 68a) for several nuclei in the 
2s-ld shell. 


10 1 ^- 

calculations for Be and '^0 have been carried 
out for only one value of b=2.09I'- He calculated binding 
energies per particle for ^^Be and are -0.77 MeY and 

-3.33 MeY compared to the experimental values of -6,5 MeY 
and -7.5 MeY respectively. Thus, "^Be is considerably under- 
bound under the s-ssiuiption of a spherically syiometric HF field. 
Our interest in doing the HF calculations for these nuclei is 
not only to study their ground state static properties but 
also to employ the HI' single pa.rticle energies and wave 
functions obtained in these calculations in a configuration 
mixing calcula.tion described in a la.ter Oliapter. In view pf 
this, viQ shall only present the results of single particle 
energies and expansion coefficients fcr these nuclei. Ihe 
neutron single particle HF levels are given in fable IY.5. 

In the case of "^'Ee » lPyy'2 lies below the lP^y'2 1®^®^ 

in absence of centre-of-mass correction, an effect similar 
to that encountered in 0. Figure IV. 2 shovrs the variation 
of neutron single particle HF levels with the nucleon number 
A. Proton energies would show similar variation. As seen 
from the figure, the levels become more and more bound as 
the nucleon number increases and the effect is largest for 
the lov/est ls^y '2 1®'^®^* 
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Table IY.3 

HF single particle energies of various nuclei for 
effective Yale interaction obtained without 
coulomb and centre-of-raass corrections. 


Nucleus 

States 

10^ 

364. 

(Me 7) 

"Ob 

(Me 7) 

16 q 

^8 
(Me 7) 

1=1/2 

-26.99 

-47.10 

-44.85 

l*' 3/2 

- 4.91 

- 20.99 

-18.95 

151/2 

- 2.56 

-16.94 

-13.78 

115/2 

- 1.38 

- 0.66 

0.67 

2=1/2 

2.50 

0.43 

1.59 

II3/2 

- 

3.10 

4.70 


SINGLE PARTICLE ENERGY (MeV) 
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The single particle energies can also "be obtained 
directly by confuting the self-consistent total energy’' E 

ex 

of the A+1 (for unoccupied states) or A-1 (for occupied states) 
system and tahing the difference of and (the ground 
state energy of the nucleus wath A nucleons) ♦ Such a calcul- 
ation has been C;arried out for ^*^0 • v;here a single particle 
was removed from the states 15^/2* ^-'3/2 ^^1/2 

to 2 s ^^2 state in turn and it is assiimed that resulting states 
15 17 

of '^Q and 0 still remain spherical. Ihiis, for the occupied 

jL ^ 

levels of 0, the energy of the i level is given by 


, = S(^^O) - E^(^^O) , 


(4.1) 


"tjix 

where i indicates a hole in the i level. Ihe unoccupied 
levels are similarlj’’ given by 


e. = E^(^^O) - E(^^O) . 


(4.2) 


The difference between the single particle energies obtained 
in this way and the corresponding energies obtained in the 
HE calculation is the rearr-angement energy of the system. 

The results are presented in Table IV. 6, 


In a similar manner, particle-hole excited states 
of 0 are calculated assuming that they are still spherical. 
Some such states are given in Table IY.7- The first excited 
0"^ state of ^^0 is experimentally observed at 6.05 MeV and 
is expected to arise from two particle-tTO hole (2p-2h), 
four particle-four hole (4p-4h) excitations whereas the 
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TalDle IY.6 


Comparison of single particle neutron energies 
obtained directly ?d.tli eigenvalues of -^^O for 
effective Yale interaction, 13=2.091'. 

= + (SQ(|;j) - Eq(GS)), the plus 'being for 

particle states and mini’s for hole states* 


State 

Total 
energy of 
the state 
E^CMeV) 

Excitation 

energy 

Bex 

Eigen- 

value 

Rearrange- 
ment energy 

*^al3"'®ex i 
(MeV) 

G-.S* 

-57.835 

0.000 



(lSl/2) 

-17.061 

-40.773 

-45 . 596 

4.822 

(lP 3 / 2 )”^ 

-38.734 

-19.101 

-21.208 

2.107 

(lPx/2^"^ 

-43.572 

-14.263 

-15.755 

1.490 

( 2 ^2 ) 

-57.921 

- 0.086 

- 0.528 

0.442 
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Table IV. 7 

T 

Particle-hole states in 0 lor effective Yale inter- 
act! on* b=2.09Y. S is determined from the 

Jx 

difference of total energies of excited 
state and ground state. 


St a- 

fee 

Total energy 
Eo{r:c.v) 

Excitation 
ene rgy 

Eg^/MeV) 

Ground State 

-57.S35 

0.000 

(Ip^/g)"^ 

(1^5/2) 

-50.436 

7.399 

(Ip^/g) 

( 2 S 1 / 2 ) 

-43.380 

13.455 

(lPl/ 2 ) 

(1^5/2)^ 

-42.100 

15.655 

(ip^/gr^ 

{2 Si/2)2 

-30.651 

27.184 


( 2p^ ^2 ) 

13.669 

71.504 


^^^1/2^ 

13.894 

71.729 
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lowest calculated 0"^ state comes at 15.65 MeV. Thus, it is 
necessary to introduce departure from spherical symmetry 
or residual interaction to describe this state. In fact, 

it is an ezperimental fact that all nuclei which are spherical 

1 fi 

in ground state, such as 0, exhibit deformation in their 

excited states and vice versa. Ihus, the excited states of 

16 

0 would not be described well under the assumption of a 

zero order spherical field. Kelson (Ken 65) and Bassichis 

and Eipka (Bas 65) have satisfactorily reproduced the excited 
16 

states of 0 by carrying out deformed HP calculations for 
the excited particle hole states and using projection of 
angular momentum techniQ[ue. They have obtained lowest 0"^ 
state very near the experimental value. 

IV. 2b Hartree-Poclc Wave Punctions and Density 

After a, sufficient number of iterations have been 
performed the single particle wave functions are determ- 

ined frora a final set of coefficients c“ by substituting 
them into eq.uation 3.20. These expansion coefficients are 
given in Tables (IV. 8 - IV. 13) for the nuclei %e, ®Be, 

^^0, ^^Be and ^*^0 respectively. The results are presented 
for the values of b v/hich give maximum binding (b=1.76P for 

^He, b=2.09P for all other nuclei). We shall discuss here 

1 1 ^ 

the nature of the wave functions for 0 alone and only make 
remarks on the structure of wave functions for other nuclei. 
The radial HP v/ave function (see equation 3.19 for 
detailed expression) for the P 3/2 ^1/2 



TalJle IV. 8 

Expansion coefficients for "o = 1,76F 

olDtained^ with effective Yale interaction. 


tes 


neutron 



proton 


a \ 

1 

2 

3 

1 

2 

3 

1 

0.9684 

0.2224 

0.1129 

-0.9694 

-0.2185 

-0.1123 

2 

-0.1573 ^ 

0.8959 

-0.4156 

0.1527 

-0.8940 

0..4213 


3 0.1956 -0.3847 -0.9025 0.1924 -0.3912 -0.8999 

1 0.8055 -0.4828 0.3435 0.7939 -0.4966 0.3510 

. 2 0.5920 0..6800 -0.4325 0.6078 0.6658 -0.4327 

* 

3 0.0247 -0..5518 -0.8336 0.0188 -0.5569 -0.8304 

1 0.7137 -0.5986 0.3636 0.7042 -0.6060 0.3700 

r 2 -0.6931 -0.5286 0.4901 -0.7020 -0.5164 0.4905 

3 0.1012 0.6018 ' 0.7922 0.1062 0.6051 0.7890 


Tab le IV, 9 

Expansion coefficients CJ^( ) for ®Be , , b = 2.09E 
obtained with effectiTe Yale interaction. 


state 

JS 


neutron 


proton 


i ^ 


1 

2 3 

1 

2 

3 


1 

0.9191 

0.3255 0.2220 

0.9204 

0.3220 

0.2219 

0 1 

2 

-0.2294 

0.9003 -0.3701 

-0.2216 

0.8970 

-0.3825 


3 

0.3204 

-0.2893 -C.9021 

0.5222 

-0 . 30 28 

-0.8969 


1 

0 . 89 30 

-0.2663 0.3629 

0.G740 

-0.3059 

0.3775 

1 1 

2 

0.3953 

0.8494 -0.3497 

0 . 4401 

0.8277 

-0.3482 


3 

0.2151 

-0.4557 -0.8637 

0.2060 

-0.4705 

-0.8580 


1 0.8025 -0.4665 0.3721 0.7816 -0.4898 0.3862 

1 I 2 0.5921 0.6990 0.3997 0.6203 0.6750 -0.3994 

3 0.0739 -0.5410 -0.8373 0.0651 -0.5518 -0.8314 
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Tablo IV. 1 0 

Expansion ooef ficients C^(2, j ) for b = 2.09E 

obtained with, effective Yale interaction. 



states 


neutron 


proton 


i 

3 

aV 

1 

2 3 

1 

2 

3 



1 

0.9357 

0.3140 C.1608 

0.9379 

0.3083 

0.1594 

0 

1 

2 

2 

-0.2430 

0.9041 -0.3516 

-0.2294 

0.8952 - 

■0.3820 



3 

0.2557 

-0.2899 -0.9223 

0.2604 

-0.3217 - 

■0.9103 



1 

0.9672 

0.1215 0.2233 

0.9685 

0.0844 

0.2344 

1 

3 

2 

2 

-0.0115 

0.8983 -0.4392 

0.0317 

0.8915 - 

■0.4518 



3 

0.2539 

-0.4222 -0.8702 

0.2471 

-0.4450 - 

■0.8608 


1 

2 


1 


1 

2 


3 


0.9667 0.0354 0.2414 0.9655 0.0453 0.2563 

0,0296 0.8991 -0.4367 0.0768 0.8912 -0.4471 

0.2543 -0.4293 -O.G666 0.2487 -0.4513 -0.8570 
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'Salle 17 . 11 (a) 

X 6 

Expansion coefficients 0^(t j ) for ' Og, 1 = 2.09E 
obtained ^'/ith effective Yale interaction. 


sta 

,tes 


neutron 


proton 

1 i 

\n 
® \ 

1 

2 3 

1 

2 3 

° 1 

1 

2 

3 

0.9472 

-0.2532 

0.1966 

0.3033 0.105S 

0.9065 -0.3330 

-0.2:939 -0.9354 

0.9497 

-0.2361 

0.2056 

0.29 6 5 0.1006 

0.8892 -0.3919 

-0.3484 -0.9145 


1 

0.9282 

0.3214 0.1876 

0.9341 

0.3034 0.1880 

i-i 

2 

-0.2118 

0.8706 -0.4441 

-0.1857 

0.8629 -0.4700 


3 

0.3060 

-0.3725 -0.8761 

0,3048 

-0.4042 -0.8624 



1 

0.9491 

0.2434 0.1998 

0.9555 0.2102 0.2068 

1 1 

2 

-0.1294 

0.8799 -0.4573 

-0.0889 0.8741 -0.4775 


3 

0.2872 

-0.4081 -0.8666 

0.2811 -0.4379 -0.8540 



1 

0.9486 

0.0570 

0.3143 

0.9460 

0.0082 

0.3240 

2 1 

2 

2 

0.0931 

0.0165 

-0.3890 

0.1242 

0.9142 

-0.3858 

3 

0.3024 

-0.3933 -0.3660 

i 0.2994 -0.4052 

-0.8638 




Sable 17.11(1:)) 


.Expansion ooefficients ) for ^^Og, b = 2.09E 

obtained in the KP oalo-Ll ati on mth effective 
Yale interaction without coulomb and 
centre-of-mass corrections. 


'i 

U' 

3 


1 

2 

3 



1 

0.9475 

0.3047 

0.0969 

0 

1 

2 

2 

-0.3008 

0.9522 

-0.0529 



3 

0.1084 

-0.0210 

-0.9939 



1 

0.9414 

0.3C32 

0.1476 

1 

3 

2 

2 

-0.2547 

0 . 9 261 

-0.2782 



3 

0 . 2210 

-0.2243 

-0.9491 



1 

0.9570 

0 . 2498 

0.1476 

1 

1 

2 

2 

-0.1953 

0.9309 

-0.3088 



3 

0 . 2146 

-0.2666 

-0.9396 



1 

0.9789 

0.0328 

0.2017 

2 

5 

2 

2 

0.0407 

0.9360 

-0,3496 



3' 

0.2002 

-0.3504 

-0.9150 



1 

0.9390 

-0 . 2555 

0.2301 

2 

3 

2 

2 

0.3302 

0.8569 

-0.3958 



3 

0.0960 

-0.4476 

-0.8891 
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le Me IV. 1 2 

Expansion coefficierits for b = 2.09 E 

obtained with, effective Yale inter.&otion* 

states neutron proton 



;] 

a 

< 

2 

3 

1 

2 

3 



1 

0.9339 

0.3023 

0.1909 

0.9363 

-0.3044 

0.1751 

0 

1 

2 

2 

-0.2150 

0.9014 

-0.3758 

-0.2277 

0.9059 

0.3570 



3 

0.2857 

-0.3099 

-0.9063 

0.2673 

-0.2944 

-0.9175 



1 

0.9548 

-0.0711 

0.2837 

0.9601 

-0.0621 

0.2725 

1 

3 

■2 

2 

O.I9O8 

0.8911 

-0. 4118 

0.1760 

0.8917 

“0.4169 

*■ 


3 

0. 2280 

-0.4483 

- 0.3643 

0.2171 

-0.4483 

-0.8671 


0.9133 -0.2546 0.3178 0.9472 -0.1594 0.2780 

0.3759 0.8272 -0.4177 -0.2708 0.8621 -0.4283 
0.1566 -0.5010 -0.0512 0.1714 -0.4810 -0.8598 

1 0.7076 -0.5740 0.4121 -0.7288 -0.5548 0.4013 

2 I 2 -0.7066-0.5660 O.4247 0.6847 0.5988 -0.4155 

3 0.0105 0.5917 0.8061 0.0097 -0.5776 -0.8163 



109 


Table IV.13(a) 

Expansion coefficients for b = 2.09? 

obtained with, effective Yale interaction. 



states 


neutron 


proton 

1 



\ ^ 

2 3 

1 

2 3 



1 

0.9420 

0.3185 0.1062 

0.9482 

0.3035 0.0942 

0 

1 

2 

2. 

--O.2715 

O.9O88 -0.3169 

-0.2525 

0.8996 -0-3564 



3 

0.1975 

-0.2697 -0.9425 

0.1929 

-0.3124 -0.9296 



1 

0-9133 

0.3552 C.I99O 

0.9199 

0.3446 0.1872 

1 

3 

2 

2 

-0.2368 

0.8610 -0,4502 

-0.2256 

0.855'5 -0.4662 



3 

0.3313 

-0.3640 -0.8705 

0.3208 

-0,3866 -0.8647 



1 

0.9268 

0.3091 0.2133 

0.9361 

0.29 16 0.19 6 6 

1 

1 

2 

-0.1850 

0.8700 -C.457I 

-0.1693 

0.8637 -0.4748 



3 

0.3269 

-0.3842 -0.3635 

0.3083 

-0.4112 -0.8579 



1 

0.9392 

0.1429 0.3121 

0.9411 

0.1614 0.2972 

2 

5 

■5 

2 

-0.0052 

0.9151 -0.4030 

-0.0277 

0.9127 -0.4078 



3 

0.3432 

-0.3770 -0.3603 

0.3371 

-0.3755 -0..8633 


Table 


Expansion coefficients 
obtained in the H? calc 
Yale interaction v 
centre-of-mac 



JL J-\J 



2 3 

0.2991 
0.9541 
0.0126 

0.3083 
0.9287 
- 0.2026 


0.2701 

0.1442 

0.9324 

- 0.2856 

- 0.2400 

- 0.9474 

0.0727 

0.1894 

0.9382 

- 0.3460 

- 0.3383 

- 0.9189 

- 0.1589 

0.2050 

0.8931 

- 0.3872 

- 0.4209 

- 0.8989 


0.0903 

- 0.0152 

- 0.9958 

0.1429 

- 0.2593 

- 0.9552 


states are plotted in Figures (IV. 3 - IV. 5). Ihe proton v;ave 
functions are only slightly different from these as can he 
noticed from the erpansion coefficients for the neutron and 
protons. Pure harmonic oscillator wave functions have also 
been calculated for a value of the parameter h which gives 
the same r.m.s. radius as obtained from HF calculations> and‘tties& 

are plotted in Figures (IV. 3 — IV. 5). From these figures we 

"1 ^ 

notice that the occupied orbitals for °0 can be well approx- 
imated by the hnrnonic oscillahor functions ?/hile the struct- 
ure of the unoccupied orbitals is quite different. The un- 
occupied orbitals are damped inside and have large amplitudes 
at large distances. This is true for all the three states. 
Similar effects are observed for the wave functions of other 
nuclei. The reason for a value of b different for HF and 
harmonic oscillator wave functions for the same value of the 
r.m.s. radius is obvious from Tables (IV. 8 - IV. 13). In other 
words, strong mixing of various orbitals is responsible for 
this difference a,s well as the nature of the unoccupied 
orbitals. From Figures (IV. 3 “ vfe also see that there 

is an overall shift of the peals of the EF wave functions with 
respect to the harmonic oscillator ones to larger values of 
r. In the spherical HF calculations one mixes different 
orbitals corresponding to the far-lying states (n mixing) 

These states have vreaker binding and are associated with 
larger separation from the centre of potential- The range 

in which nucleons move is considerably large compared to 
pure states. Therefore, the average separation of the 








F!G.nz:-5 



nucleins froin the centre of the potentisil is increased, result- 
ing in the shift of peaks. For ^^0, using the HF values of r.m.s, 
radius the calculated value of h for the oscillator functions 
turns out to be 1.699F. However, for this value of b the two 
wave functions do not overlap exactly. Best overlap is obtained 
for the tabulated b values (lable IY.14) , For ^e nucleus, the 
calculated value of b is 2.424F. This large b indicates weak 
binding as can be verified from Table IV. 1, 

HF single particle proton orbitals with and without the 
correction terms have been plotted in Figures (IV. 6 - IV. 8) for 
some of the single particle states in 0 and ^He. As seen from 
the Figures (IV, 6 IV. 7), the correction terms have very little 

effect on occupied orbitals in 0 but the ■unoccupied orbitals 
have been significantly affected and the centre-of-mass correction 
is much more important than the coulomb correction. Figure IV. 8 
shows the proton orbital in ^He with and without the centre-of-. 
mass correction. 

.4 

HF, matter and charge densities for light nuclei ( He, 

®Be, ^^0) are shown in Figures (IV. 9 - IV. 10). Both the 

16 

matter and charge densities of 0 show a dip at the centre. 

Such a dip is predicted from electron scattering experiments. 

IV.2c Hartree— Fock Potential 

The matrix elements <ntUjn^ of the HF single particle 

16 

potential U are given in Table IV. 15 for 0 neutron states for 
b=2.09F . Since < nj U |n *) = < n'^ j U {n > , only n ^ n* values are 
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fable IV. 14 


Values of the oscillator parameter ^b) giving 

16 

the best overlap of the EF orbital. for ©0 
obtained with effective Yale interaction and 
harmonic oscillator wave functions. 


Xk-bital 

bCF) 

lSi/2 

1.612 

IP 5/2 

1.625 

^^1/2 

1.667 







E? ■ 7 • p state EP piroton waTe ftoKjtions for ^^0 obtained witb 
®ff©otiTe Yale interaotioa. Other details are saiae a® In 

’ kc. I¥.S. 
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ga ble IY.13 


Matrix 

for 


elements ^xi jtJ(Pi)|n*^ of the HE potential 
0 , b = 2 . 09.2 obtained with effective Yale 
interaction for neutron states. 



n^ ^ 

1 

2 

3 


1 

-25.5188 



0 1/2 

2 

-13.6123 

13.3058 



3 

- 6.4-951 

1.6423 

25.7619 


1 

- C.1265 



1 3/2 

2 

- 7.5401 

15.6957 



3 

- 7.1378 

3.2100 

25.5869 


1 

3.9921 



1 1/2 

2 

- 5.2082 

17.8884 



3 

- 6.5334 

4.4390 

26.8043 


1 

IO.C 8 I 5 



2 5/2 

2 

- 1.8058 

16.9826 



3 

- 5.4234 

4.4253 

24.0448 



given. The potentials ) iiavo also been evaluated using 

expression (3.62) for different single particle neutron states 
(ig) of ^^0 for b=2.09P. These potentials have been plotted in 
Figures (lY.ll - IV .13) for the s^y'2^ ^3/2 ^1/2 

states respectively in 0, The fact thef HP potential is non- 
local and sta.te-dependent is obvious from the plots. A local 
potential would give a delta function at r = ri Potentials for 


and p. 


states are veiy similar. The HP single particle 


potentials for other nuclei are expected to show similar behaviour. 


The sta.tic limit of the HP potential has been calculated 
using equation (3.65). As seen from Figure IV. 14, this potent- 
ial does not resemble any of the standard local potential forms 
such as harmonic oscillator, Woods-Saxon, etc . The static limit 
of the potential is most attractive at about r = 1.8F and becomes 
repulsive beyond r = 3.6F. It is not clear what meaning this 
local potential can have. 

We have also attempted to find out the parameters of 
a local potential which would give the same single particle 
energies as obtained in the HP calculation. The local potent- 
ial employed is of Woods-Saxon form with spin-orbit and coulomb 

terms added to it. The detailed shape of the potential is given 

16 ‘ 

in Appendix A. Calculations have been carried out for 0. The | 
range parameters r^^ and r^ of the potential were kept fixed at i 
a value of 1.25F, approp mate for ^^0 and the strengths V^j and , 
V of the central and spin-orbit parts of the potential are 


□ 


6-0 


8-0 


4-0 


r(F)' 


16 

iff single parijicle potential for 0 


;ive Yale interaction. 


obtained 


fl J V^1VI ^¥ J 



FIG.EI. i 2 P5/2 ***'®^*® ® siaigl® pairUol® pottaiial for 0 obtalaat 
wiiai «ffec1:i¥» Xalo iatowustion* 




Vstdt 


0-0 


■ 0*2 



r(F) 

FiG. E- 14 Static limit of tlie non-local HP p< 
witli offoctiv© Tala interaction# 
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varied so as to reproduce the HP single particle levels. Por 
the Ip^^gC 21.21 HeY) and lp^^ 2 (l 5.75 MeY) neutron levels of ^^0, 
the strengths are fotmd to he Y„ = 53.206 MeY and Y = 8.834 

ijt SO 

MeY. The corresponding single particle wave functions are 
plotted in Figures (IY.15 - IY.16) and also compared with the 
HP wave functions. It is seen that the tyro wave functions (HP 
and equivalent local potential) agree pretty closely in the 
region r43»5P>htit the HP wave function has a iiuoh longer tail. 
Hon-local correction to the wave function (see Appendix A) has 
also been applied. The lP^y'2 functions for different 

values of the non-locality range parameter a^^ are shown in 
Figure IY.16. As the parameter a^^ is increased the tail of the 
wave function approaches the HP wave function, hut the disagre- ^ 
ement in the inner region increases. 

It should he mentioned that the two— hody coulomb and 
centre-of-mass matrix elements, used in the calculations reported 
in this section, have been calculated with a maximum relative | 

1,= 2. This has caused some error in our calculations since [ 

coulomb force is a long-range force and there would be signif— i 

icant contributions coming larger i values, later on, 

■■ . f 

coulomb and centre-of-mass matrix elements were calculated with 

a maximum relative i = 10 (maximum allowed I, in the smaller I 

configuration space) and some of the HF calculations for effect—! 
ive Yale interaction were repeated with Idle se new coulomb and 
centre-of-mass matrix elements, A few revised HF results have | 
been added to Appendix B, It is seen that static properties of | 







m 

o 

# 


m 

1 




m 

3 

o 


§5g 

© O -r* 

S m 4 » 

IN Q 

©NS 

I S « 

N t! fe 

tab P 


3 03 'tS 
^ »d § 

« i PN 


^ § a 
I :^ © 

m |£| 

H H ^ 
-H 0 « 

# 0 D 

-P « 

O P# ® 

1 ::% ' S 

■3 5 8 
g « . 

0 fc 

.rig 

^ «r 4 P 

. m 0 

%| 4 » 

s ,0 ^ 

0 m 

“ • 

'a m m 


u O' 

© © 

■ «rl 


4* 

m 

# 

m ^ 

0 

4» • 

3 

« O" 

'Cl 



0 

0 0 
H H 

3 

g 

3 g 

1 I 

g. g 

“ a 

N o 
o o 

11 

tf 0 
4* r4 


m %4 
# 

> • 

5 I 


i» ca 

s \® 

P HN 

^ ^ i 

i 0 ' 4 |fi| 

■** 5 •»• 

©a © 

P 

C>i CM H 
\ \© 
m lA-i* 
Pi 

1-4 f-4 N 


■P © 

® P 
m 44 

S >© 

Pi © 

s a 

m 

• 4^ 

§ m M 
0 

m 

O -H 


m 

rl 


0 


1^ ■ 


Tii. 



ti' 

0 

p 

# 



id ' 

m'l 

# 

1 

*3 

0 1 

$w| 

3 

pS|l 

4i 

OI 


*; 

, 0 

m\ 






131 


nuclei are not -very much affected. Ihese new coulomh and centre- 
of-mass matrix elements have been used in the HP calculations 
with Sussex interaction reported in next section. 

IV. 3 Sussex Interaction Results 

HP calculations with Sussex relative matrix elements : 
have been done for the nuclei "^He, , ^^0 and ^®Ga. Two- 

body matrix elements occurring in the HP equations are calculated 
using tables of relative matrix elements given in reference 
(Elt 68a). The tables given in reference (Elt 68a} are not 
complete since only the matrix elements between the states 
whose oscillator energies [i.e. ( 2n+E— differ by not more 
than 2h3UJ are tabulated. Some other matrix elements, such as 
off-diagonal matrix elements of the type <^n''|V|n y when the 
radial quantum numbers differ by more than one, i.e. (n -n) ^2, ; 
are also needed in the present work and have been calculated 
using a simple recurrence formula (Elt 68a) from the tables ; 

given by Elliott et. al. (see Appendix G). The complete set ; 

of matrix elements used in our calculation i n given in Tables I 
IV. 16 for some of the relative states. The specific form of | 

the two-body Sussex interaction is not known but is not needed | 
in any part of the calculation. 

Resiats of HP calculations obtained with Sussex inter- : 
action are tabulated in Tables (IV. 17 — IV, 19) for the nuclei 
^He, and ^^0 and for several values of b. We have not 





T able IY.16(a) 


Sussex relative matrix elements for channel, 
h = 2. OF, Only rJ ^ n values are given since 
<n ^3^1 Yjn^S 3 _>=<n 5 s^j Yjn^ 


n 

n' 1 

2 

3 

4 

5 6 7 

1 

-5.58 





2 

- 5.22 

-5.11 




3 

-4.41 

-4.24 

-3.76 



4 

-3.45 

-3.35 

- 2.93 

-2.53 


5 

-2.47 

-2.46 

- 2 . 21 

-1.89 

-1.61 

6 

- 1.56 

- 1.62 

-1.53 

-1.34 

- 1.13 -0.89 

7 

-0.74 

-0.87 

-0.89 

-0.83 

-0.67 -0.43 -0.07 




Table 

IV.l 6 (h) 


Susse^f relative 

matrix 

3 

elements for channel, 


h = 2. OF, Only 114 n 

values 

are given since 




= <'n^P 

3^|V|n'\>. 

n 

n' 1 

2 

5 

4 

5 6 7 

1 

1.20 





2 

1.34 

1.85 




3 

1.40 

1.97 

2.48 



4 

1.41 

2.13 

2.70 

3.14 


5 

1.42 

2.27 

2.87 

3.26 

3.39 

6 

1.47 

2.36 

2.95 

3.20 

3.22 3.21 

7 

1.53 

2.4 0 

2.38 

3.00 

3.04 0.00 0.00 


17.16(c) 


Sussex relative matrix elements for ciiannel, 
"b = 2. OP. Only n'*Cn values are given since 
<n' ^P^ j 7 fn^P^) = (n^p^ | 7 jn' 


a 

1 

2 

•z 

J 

4 

5 

6 7 

1 

0.86 






2 

1.21 

1.78 





3 

1.46 

2.17 

2.71 




4 

1.66 

2.48 

3.09 

3.54 



5 

1.81 

2.72 

3.37 

3.84 

4.14 


6 

1.93 

2.88 

3. 56 

4.01 

4.30 

4.50 

7 

2.02 

2.99 

3.65 

4.08 

4.39 

4.63 0.00 




Sable 

17.16(d) 



Sussex 

relative 

matrix 

/ 

elements for 

5 

Dg obannel. 


Td = 2 

.OP. Only n4 n 

. 'Z 'Z. 

values 

are given since 




n 

xJ 1 

2 

3 

4 

5 

6 7 

1 

-1.63 






2 

-1.43 

-2.01 





3 

-1.19 

-2.11 

-2.46 




4 

-1.09 

-1.92 

-2.33 

I 

i 

a\ 

o 



5 

-1.04 

-1.59 

-2.15 

-2.44 

-2.74 


6 

-0.93 

-1.34 

-1.87 

-2.34 

-2.62 

-2.86 

7 

-0.77 

-1.16 

-1.6 2- 

-2.15 

-2.51 

-0.00 -0.00 


1 
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Table IY.17 


Oaloulated HF properties for obtained mth Sussex 

interaction. 


b(F) 

1.5 

1.8 

2.0 

2.2 

Experi- 







mental 

B.E./a (Me¥) 

- 3.602 

- 3.263 

- 3.134 

- 2.947 

- 7.07^ 

Matter 

r.m. s* 

. 1.695 

1.827 

1.898 

1.999 


radius 

(p) 




Charge 

r*m* s. 

1.699 

1.831 

1.901 

2.002 

1.67^ 

radius ' (F) 





Single particle 
energies (MeV) 






Is^ 

2 

neutron 

-24.05 

-21.24 

- 20.15 

-18.75 


proton 

-23.21 

-20.47 

- 19 . 41 

-18.04 

-20.4+0.3° 


neutron 

6.38 

4.62 

3.94 

3.42 

0.96'^ 

IP 3 

2 

proton 

7.36. 

5. 46 

4.72 

4.12 

1.97'^ 


neutron 

8.86 

6.05 

4.99 

4.18 

4.36'^ 

IPl 

2 

proton 

9.76 

6.83 

5.70 

4.81 

5.37'^ 

2s^ 

neutron 

15.53 

8.65 

6.83 

5.57 


proton 

14.47 

9.40 

7.50 

6.18 



neutron 

2.49 

1.43 

1.05 

0.76 


IP 3 -IP 

1 













splitting 


1.36 

0.98 



0.68 


(MeV) 

proton 

2.41 



a. lattaucii et, al. (Mab. 65) 

b. Hofstadter et. al. (Hor 67) 
G. Eiou (Eiu 65) 

d. Vashakidze et. al-. (Vae 68) 
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Tab le IV. 18 


Oaloulated HF properties for obtained with Sussex 

interaction . 


b(p) 

1.5 

1.8 

2.0 

2.2 

Experi- 

mental 

B.E./a (MeV) 

- 2.309 

- 2.107 

- 1.992 

- 1.828 

- 7.68^ 

Matter r .m. s. 
radius (F) 

2.227 

2.404 

2.489 

2.609 


Okarge r.m.s. 
radius (F) 

2.238 

2.419 

2.504 

2.625 

2.47^ 

Single particle 
energies (MeV) 






neutron 

-47.73 

-42.03 

-39.67 

-36 , 51 


Is^ 

-2 proton 

-44.67 

-39.25 

-36.99 

- 33.94 

-34.2+2° 

-33.0+3.5 

neutron 

IP 3 

proton 

- 15.34 

- 13 . 49 

-12.79 

- 11.90 


-12.54 

- 10.92 

- 10.30 

- 9.52 


neutron 

-11.98 

-11.78 

-11.77 

-11.47 


IPl 






proton 

- 8.84 

- S.85 

- 8.91 

- 9.25 


neutron 

5.04 

2.76 

2.34 

2.18 


2s^ 

proton 

7.71 

5.00 

4.49 

4,17 


T T neutron 

IP 3 -IPP 

3.36 

1.72 

1.02 

0.43 


2 2 






splitting 
(MeV). proton 

3.70 

2.07 

1.39 

0.29 



ja. Mattauoh et. al» (Mah. 65) 
b. Hofstadter et. al. (Hor 67) 
G. Biou (Eiu 65) 
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gclle IV. 19 


1 

Calculated HF properties for Og otitained with Sussex 

interaction. 


b(F) 

1.5 

1.8 

2.0 

2.2 

Experi- 







mental 

B.E./a (MeV) 

- 5.111 

- 4.910 

- 4.850 

- 4.772 

- 7.98®- 

Matter r. 

m. s . 

2.217 

2.325 

2.380 

2.478 


radius r^ 

,(B) 



Charge r. 

m, s. 

2.227 

2.337 

2.391 

2.479 

2.71^ 

radius r^ 

.w 



r^.X B.E. 
(MeY-E) 

/A 

11.231 

11.416 

11.543 

11.825 

^ 11° 

Single particle 






energies 

(MeV) 







neutron 

-57.16 

-55.05 

-51.53 

-48.45 

tA -51 

ISi 

2 

proton 

-53.04 

-49.16 

-47.51 

-44.65 

-34+3.5® 

-44+2 

IP 3 

neutron 

-28.82 

-26.57 

-25.77 

-24.44 

— 21. 81*^ 

-18+2.5® 

proton 

-24.86 

-22.82 

-2 2.11 

-21.17 

2 




-19+1 

ipi 

2 

neutron 

-19.90 

-19 .18 

-19 . 20 

-18. 90 

-15.65*^ 

proton 

-16.06 

-15.55 

-15.61 

-16.44 

-15+2® 





-12.4+1 

2 Si 

2 

neutron 

2.12 

0.15 

0.10 

0.31 

- 4.15^ 

proton 

5.76 

3.38 

3.15 

3.12 


Ip^-lPl 

neutron 

8.92 

7.39 

6.57 

5.54 

6.16^ 

2 







splitting 
(MeV) proton 

8.80 

7.29 

6.50 

4.73 



a. Mattauch et. al. (Mah 65) 
h. Hofstadt*!' et, al. (Hor67) 
o. Mgeli et. al. (Ani^69) 

d. Cohen et. al. (Con 63a, 631) 

e. Eiou (Eiu 65) 
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presented the results for h = 1.4F and h = 1.6F, since they are 

close to h = 1.5P results. Binding energy per particle vs. b 

curves (Figure IY.17) are very nearly straight lines for all 

these nuclei (^He, and ^^0). ¥e find that the calculated 

binding energies for 0 is about (30 — 40 percent) less than 

the corresponding ezperiinental value and still poorer agreement 

with the experimental data has been obtained in the case of "^He 
12 

and 0. As described in Section IY.2a, binding energy result 
12 

for 0 will improve if the constraint of spherical symmetry 
is relaxed. Moreover, since minimas in the binding energy vs. 
b curves could not be obtained till b = 1.5F, the calculations 
should have been extended to lower values of b in order to loiow 
how much maximtua binding one can obtain in such a calculation. 
The matrix elements in reference (Bit 68e) have been calculated 
only up to b = 1,4F since lack of experimental phase shift data 
above energies )> 300 MeY imposes a lower limit on the value of 
b up to which matrix elements can be determined accurately (see 
Section II. 3) . However, as pointed out by Elliott et. al, the 
matrix elements for lower b values can be obtained by a 
polynomial , fit to the matrix element vs. b curve in the known 
region and extrapolation of the fitted polynomial to smaller. b 
values, but the reliability of the matrix elements so obtained 
can be questioned. Moreover, even for the smaller configuration 
space included in the present calculation and for b^l.4F, the 
relative matrix elements of Y involving energies E> 300 MeV will 
also enter. The second and higher-order terms (see expansion 
2 . 56 ) for such matrix elements are quite important since the 



' binding ene; 
,ssex interao' 
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a'uxiliar;r potential can be chosen to approximate the true potent 
ial only below 300 IieV. In view of this, the calculations have 
not been carried out for smaller values of b. 

As seen froti I'able IV..19 , the single particle energy 
levels in 0 are overbound compared to experimental data. The 
coulomb correction shifts the proton levels relative to neutron 
levels but the relative separation between any two proton levels 
is very nearly the same as the corresponding separation between 
the neutron levels. Angeli et. al. (Ani 68), on the basis of 
experimentally measured cross-sections, have suggested the fol- 
lowing relationship betv/een the binding energ^^ per particle and 
matter r.m.s, radius for ^^0, 

B.E./A (MeV) X 11. (4.5) 

Our HP calculations for ^^0 agree with this. 

Assuming that the Sussex interaction V can be broken 
up into central, spin-orbit and tensor parts, the matrix element 
of any particular component can be expressed in terms of the 

elements of the full potential as described in Appendix 0 
Such a calculation has been carried out and the results are pres 
ented in Table IY.20 for some of the relative states of nucleon- 
nucleon pair. In order to test the effects of each of these 

Various ground state properties hi calculations have 
been repeated with the contributions of 
(i) only the central part. 
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gable IY.2C(a) 


Central, 

spin- orbit and 

tensor components of ■ 

Sussex 

relative 

matrix 

elements 

< n|V{n>, 

diagonal in 

osci- 

llato 

r radial quantni'i 

. number n 

for b = 2-OP. 

Radial Gliannel 

Quantem 

Niimber 

Central 

Spin-orbit Tensor 

Total 


U 

-0.13 

0.60 

-2.05 

-1.58 

n=l 

^1 

-0.13 

O 

« 

o 

1.03 

1.20 


^ 2 

-0.13 

-0.30 

-0.21 

-0.64 


3p 

^0 

-0.15 

1.34 

-2.67 

-1.48 

n=2 

^P 

^1 

-0.15 

0.67 

1.33 

1.85 


d. 

-0.15 

-0.67 

-0.27 

-1.09 


^P 

•^0 

-0,03 

2.11 

-2.92 

-0.84 

n=3 


-0.03 

1,06 

1.46 

2.49 



-0.03 

1 

• 

o 

a\ 

-0.29 

-1.38 



-0.39 

0.10 

1 . 27 

0.98 

11=1 


% 

o 

1 

0 . 03 

-1.27 

-1.63 


^P3 

-0.39 

-0.07 

0.36 

-0.10 


Contd. . 
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Tabl e I V, 20(a) ( 0 o n t d . ) 


Radi al 
Quantum 
Number 

Channel 

Central 

Spin- orbit 

Tensor 

Total 



-0.46 

0.28 

1.65 

1.47 

n='2 


-0.46 

0.09 

-1.65 

-2.02 



-0 . 46 

-0.19 

0.47 

-0.18 



-0.56 

0.43 

2.05 

1.92 

11 


~C . 56 

0.14 

-2.05 

-2.47 



- 0.56 


0.29 


0.58 


-0.27 
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2 able I V. 20(1) 

Central, spin-orbit and tensor components of Sussex 
relative matrix elements <( n jV|n > , off-diagonal in 
oscillator radial iiiantum number 
for b = 2.0'f . 


Radi al 

Quantum 

Number 

Gliannel 

Central 

Spin- orbit 

Tensor 

Total 

n = 1 

^0 

3 

-0.15 

0.88 

-2.1 

-1.37 

n'= 2 


-0.15 

0.44 

1.05 

1.34 


-0.15 

-0.44 

-0.21 

-0.80 


^0 

-0.10 

1.68 

-2.46 

00 

00 

o 

1 

n = 2 ^ 



-0.10 

C.84 

1.23 

1.97 

n'= 3 


^P 

^2 

-0.10 

-0.84 

-0.25 

-1.19 


5p 

^0 

-0.13 

1.11 

-1.95 

-0.97 

n = 1 

n' = 3 

2p 

^1 

-0,13 

0.55 

0.98 

1.40 



-0.13 

-0.55 

-0.20 

-0.88 



-0.31 

0.22 

1.19 

1.10 

n = 1 

i 


-0.31 

0.07 

-1.19 

-1.43 

n = 2 


^3 

-0.31 

-0.15 

0.34 

-0.12 


Contd 
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(ii) centre;.! -rid spin-orlit part, 

(ill) centra.1 end tensor parts of tie interaction. 

We present the resnlts in ‘la’Dles (17,21 - 17.22) fox the nuclei 
^^0 and respectively .and for b = 2. OP. ihis is also the 
value of b for t.ii-cI spin-orbit split-jing for 0 is closest to 
experiment?.! v.aiue (see fable 17.19). As seen from the tables, 
major portion of t’:i.e binding energ*' conns from central part of 
the interaction v.hile the tensor force leads to saturation. In 

the Cc.se of 0 cluost all of the spiii--orbit splitting comes 

IP 

from the spiii-or’.jit part of the force. In C, tensor force 

also produces large splitting but of urong sign, fhis is cons- ; 

istent with the findings of Wong (Wog S3) and farbutton (fan 68) i 

fhe contribution to the splittings from the spin-saturated 

shoj-ls (in spin-saturo>ted shells both the members of the doublet 

j = t+ ^ are coipiletely filled), such as in °0, is positive and| 

comes mainl?’ from the spin-orbit part of the two-body potential.; 

12 

In the nuclei witl: spin unsaturated s.-ells, such as in C, | 

there is significent contribution to tl:ie splitting from the I 

tensor forces, but it always has the wrong sign. Thus, the j 

effect of only partially filling up a j = t, + i doixblet is to 
decrease the splitting, ihe spin-orbit splitting is a sensitive 

j 

function of b .and increases rapidly -with decreasing b. fhis is 
expected since the spin-orbit force dominates at high-energy . j 

40 i 

H~ p calculc,tions have also been carried out for Oa 
and also for ^*^0 in the larger configuration space which I 
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■Table IY.21 

Effect of central, spin- orbit tensor parts of 

Sussex interaction on calculated HP properties of 

for b = 2. OP. 


Only Central 


B.E./^ (MeV) 

-4.828 

Matter r.m.s. 

radius 

(P) 

2.376 

Charge r.m.s. 

radius 

(E) 

2.387 

Single particle 
energies (MeV) 

Neutron Proton 

isi 

2 

-51.11 -47.29 

IP3 

2 

-23.60 -13.96 

IPl 

2 

-23.60 - 13.96 

2s^ 

2 

0.15 3.20 

IP3-IPX 

2 2 
splitting 
(Me?) 

0.00 0.00 


No Tensor 

No Spin-orbit 

- 4.924 


- 4.917 

2.365 


2.362 

2.367 


2.365 

Neutron Proton 

Neutron Proton 

-51.78 - 

■47.82 

-51.61 -47.70 

-26.13 - 

•22.63 

-23.86 -20.40 

- 19.15 - 

■16.25 

- 23.73 -20.93 

0.14 

3.18 

0.20 3. a 

6.98 

6.38 

0..13 -0.53 
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l a}3le IV. 22 


Effect of centrals si'in-orbit and tensor parts of 
Sussex interaction on calculated HP properties of 

^^ 0 ,- for H = 2. OP. 



Only Centre!. 

No Tensor 

No Spin- orbit 

B;E,/^ (Me?) 

- 1.506 

-2.627 

-1.007 

Matter r.m.s, 

radius 

(P) 

2.594 

2.403 

2.712 

Charge r.m.s. 

radius 

(S') 

2.610 

2.415 

2.731 

Single particle 
energies (MeY) 

ileutron Proton 

Neutron Proton 

Neutron Proton 

Isi 

2 

-54.25 -31.72 

-41.49 -38.71 

-32.32 -29.90 

IP 3 

2 

-11.67 - 9.50 

-15.51 -12.93 

- 9.34 - 7.07 

IPl 

2 

-10.49 - 7.72 

- 8.74 - 6.13 

-12.54 -10.09 

2Si 

2 

2.63 4.73 

2.45 . 4.62 

2.51 4.58 

IP^-lPl 

2 2 
splitting 
(MeV) 

1.18 1.58 

6.77 6.80 

- 3.20 - 3.02 
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includes d^y '2 oscillator states of n = 1 - 3 major 

shells as v/ell. In these calculations, coulomb and oentre-of- 
mass corrections were not included in a self-consistent manner. 

Results are presented in Tables (IY.23 - IY.24) for ‘^^a and 

jL 6 

0 respectively and for b = 2. OP. In the case of the 

calculated binding energy is close to experimental value, but 
the single particle levels are somewhat overbound compared to 
experimental data. Expansion coefficients are given in Tables 
(IY .25 - IV. 28) for the nuclei "^He, ^^0 and ^^Ca respect- 

ively. Table IV. 29 shows variation of expansion coefficients 
16 

for 0 v/ith b. 


IV. 4 Discussion 

In our calculations with effective Yale interaction 
we find that binding energies for "^He and ^^0 are in fair 
agreement with the experimental data when the second-order 

term in VP is talcen into accoxint. Reasonable values for the 

8 12 

binding energies of Be and C can only be obtained if the 
HP field is permitted to be deformed. In important correction 
to the binding energ;)^ arises because of the approximate 
treatment of the energy denominator e in the evaluation of 
second-order terms in expression (2,27). The energy denom- 
inator e is defined as the difference in the energy between 
the excited intermediate state and the occupied initial state. 
In reference (Shn 67a), the intermediate states were tahen to 
be plane waves and the energy of the interacting pair was 
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T alDle IV>23 

Oalculated HF properties for obtained \»nth 

Sussex interaction vjitliout coulonb and oentre-of — 
mass corrections, b = 2. OP. 



Sussex 

Experimental 

Single particle energies 
(MeV) 

^®l/2 

-93.04 


IP 3/2 

-61.70 


^^1/2 

- 55.03 


1^5/2 

-33.17 

- 21 . 9 ® 

^®l/2 

- 25.88 

-18.2®- 

ld3/2 

-22. 20 

- 15 . 6 ^ 

B.E./A (MeT)., 

- 9-.84 


=ooul/A 

1.85 


(Me 10, 

- 0.20 


B.E./a (MeV) including c.m. 

- 8.19 

- 8.55^ 

and coul. oorreotions 



Ip^/g-lp^/p splitting (i;eV) 

6.67 


splitting (MeV) 

10.97 

7-.0° 

iL-splitting in the same major 

- 2.90 

- 1 . 1 '^ 

shell, 2s-ld (MeV) 

J.m.s. radius (P) 

2.789 

3 . 50"^ 


a. Elton et, al. (Eln 67) 

b, Mattauoli et. al- (Mah. 65) 

G. Oohen et, al. (Con 63a) 

d. Dey et- al- ( Dey-P reprint) 

e, Hofstadter et, al. (Hor 67) 
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Table IV. 24 

Single particle energies of "^^Og, b = 2. OF obtained 
in the HF oalonlation with Sussex interaction 
?hthout coulovio and centre-of-mass corrections* 


State Single particle energy 

(MeV) 


Is^ 

-51.08 

2 


Ip^ 

2 

-23.41 

IPl 

-17.13 

2 


1% 

- 1.87 

2 


2si 

- 0.68 

2 


Id^ 

3.40 


2 


Expansion ooeffioients C^(lj) for obtained 





va. th. Sussex into 

iraoti on 

j b = 1. 

SF. 



states 

neutron 



proton 



3 

\n 

a 

\ 

1 2 

3 

1 

2 

3 



1 

0.9700 0.2152 

C.112S 

0.9711 

0.2108 

0.11 21 

0 

1 

2 

2 

-0.1599 0 . 915 c - 

• 0.3704 

- 0.1547 

0-9131 - 

-0.3773 



3 

0.1829 -0.34-13 - 

■ 0.9220 

0.1818 

- 0.3490 - 

- 0.9193 



1 

0.9020 -0.3627 

0.2341 

0.8939 

-0.3785 

0.2401 

1 

3 

2 

2 . 

0.4301 0.8018 - 

■ 0.4150 

0.4472 

0.7907 - 

-0.4182 



3 

0.0371 - 0.4750 - 

-0.8792 

0.0316 

-0.4812 - 

-0.8761 



1 

0.8 3 23 -0.4349 

0,2495 

0.8239 

-0.5057 

0.2558 

1 

1 

2 

2 

-0.5502 -0.6833 

0.4300 

-0.5621 

-0.6718 

0.4823 



3 

-0.0671 0.5368 

( .8410 

0.0721 

0.5412 

0.8378 
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l alle I¥,26 

Expansion ooefficients 0^(1;} 3 for obtained 

?ith Sussex interaction, b = 2,0E, 



states 


noiitron 



proton 


h 

0 

a\ 

1 

2 

3 

1 

2 

3 



1 

0.9397 

0.3139 

0.1359 

0.9422 

0.3075 

0.1330 

0 

1 

2 

2 

-0.3004 

0.9473 - 

-0.1109 

-0 . 2888 

0.9467 - 

-0.1427 



3 

0.1635 

-0.0634 - 

-0.9845 

0.1698 

-0.0960 - 

-0.9808 



1 

0.9647 

0.2140 

0.1532 

0.9688 

0.1954 

0.1525 

1 

3 

2 

2 

-0.1582 

0.9365 - 

-0.3126 

-0.1360 

0.9335 - 

■0.3318 



3 

0.2104 

-0.2773 - 

-0.9374 

0.2072 

-0.3007 - 

-0.9309 


1 

2 


1 


1 

2 


3 


0,9675 0.193S 0.1547 0.9714 0.1801 0.1546 
-0.1451 0.9404 -0.3074 -0.1220 0.9376 -0.3258 
0.2069 -0.2750 -0.9389 0.2036 -0.2976 -0.9327 


iNol H j I rol H 




1 0.9433 0.3140 0 

2 0.3163 - 0 . 94 G 7 -0 

3 - 0.1011 - 0.0379 0 


.1079 

0.9461 

0.3069 

0.1038 

.0040 

- 0.3040 

0.0517 

- 0.0433 

.9942 

0 . 11 21 

- 0.0094 

- 0.9937 


1 0.9236 0.3447 C .1677 0.9286 0.3328 0.1639 

2 - 0.3026 0.9242 - 0.2331 - 0.2837 0.9 218 - 0.2643 

3 0,2353 - 0.1645 - 0.9579 0.2391 - 0.1989 - 0.9504 


1 0.9445 0.2882 0.1577 0.9500 0.2712 0.1549 

2 - 0.2380 0.9311 - 0.2765 - 0.2153 0.9279 - 0.3043 

3 0.2265 - 0.2236 - 0.9480 0.2262 - 0.2557 - 0.9399 
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lalD le IV. 27(b) 


Expansion coefficients Rj ) for obtained in 

Ji. O 

the HF calculation with Sussex interaction without 
coulomb and centrc--of-mass correction, b = 2. OF. 


state 

^ \ 

1 

2 

5 



1 

0.9393 

0.3224 

0.1132 


0 1 

2 

0.3275 

-0.9444 

-0.0287 



3 

-0.0976 

-0.0640 

0.9952 



1 

0.9336 

0.3216 

0.1582 


^ i 

2 

-0.2817 

0.9314 

-0.2305 



3 

0. 2215 

-0.1706 

-0.9601 



1 

0.9555 

0.2609 

0.1485 


J. 2 

2 

-0. 2136 

0.9374 

-0.2749 



3 

0.2107 

-0.2306 

-0.9500 



1 

0.9750 

0.1027 

0.1971 


2 — 

^ 2 

2 

-0.0325 

0.9428 

-0.3318 



5 

0.2199 

-0.3172 

-0.9225 



1 

0.9592 

-0.1930 

0.2065 


2 1 
^ 2 

2 

0.2618 

0.8821 

-0.3917 



3 

0.106 G 

-0.4298 

-0.8966 
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T able IV. 28 

Exp^ansion ooefficients ) for obtained in 

the HE oalciilation ’rtith Sussex interaction without 
coulomb and centre-of-mass correction, b = 2. OP. 


states 

\ 





h 2 

\ ^ 
a \ 

1 

2 

3 



1 

0.9455 

0.3139 

0.0870 


0 1 

2 

0.3211 

-0.8533 

-0.4108 



3 

-0.0547 

-0.4163 

0.9076 



1 

C.9163 

0.3781 

0.1320 


1 I 

2 

0.3995 

-0.8859 

-0, 2357 



3 

-0.0278 

-0.2687 

0.9628 



1 

0.9167 

0.3740 

0.1407 


1 1 

2 

G . 39 27 

-0.9082 

-0.1445 



3 

-0.0738 

-0.1878 

0.9794 



1 

0.9058 

0.3886 

0.1687 


5 

2 2 

2 

-0.3779 

0.9212 

-0.0926 



3 

C.1914 

-0.0201 

-0.9813 

.. 




1 

C.9316 

0.-3266 

0.1595 

2 

3 

2 

2 

-0.2938 

0.9350 

-0.1986 



3 

0.2140 

-0.1381 

-0.9670 
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a?aole l Y, 29 

1 f) 

Expansion ooefficients for “ Og obtained with Sussex 
interaction. Only j ^ ooeffioients 

are given. 


1 

d 

b(P)-~> 

ooefficient 

.si . . 

1.5 

1.8 

2.0 

2.2 



Cl 

-0 , _ -^4 

0.9729 

0.9433 

0.9162 

0 

1 

^2 

-0 .9553 

0.9703 

-0.9487 

-0.9199 



g5 

^3 

-0.9674 

-0.9936 

0.9942 

0.9855 



^1 

-0.9;J35 

0.9659 

0.9236 

0.8827 

1 

3 

C2 

^2 

-0.9243 

0.9373 

0.9242 

0.9063 



3 

^3 

-0.92.17 

-0.9544 

-0.9579 

-0.9510 




0.9983 

-0.9809 

0.9445 

0.9043 

1 

1 

2 

4 

0.90' '2 

-0.9357 

0.9311 

0-9170 



4 

-0.9075 

-0.9395 

-0.9480 

-0.9452 
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replaced Dy a coj-isacO-ilj avera.ge v'alue Zx * riiusi e was 'SiclieiL as 


e = 


+ =4) - A 


(4.4) 


?/it}aA = 20 Iio¥. Zowever, A is actually' state -dependent and 
generally larger than 20 KeV. llie approzinate tr-eatinent of 
causes an uncerts.i..it 3 '' of about 20 percent in tlie nj.triz 

elements wbicl: are the most important ones in the binding energy 
calculation. Actuallj', the energy dencnlnator should be 


-St 


hr 


AlfT-2 . v2n , _ , 

2m A A + ^2^ + ^lilh 

A 




(4.5) 


where e^.Fi’ ^2H? single particle energies 

of the initial state of the nucleon pair. This v;ould mean tioat 
the second-order terns occurring in t/ie effective interaction 
depend upon the !IiJ‘ solutions, which an turn cannot be obtained 
?dthout a pre-Knowledge of the effective anteraction. Thus, 
there is a double self-consistency proilei-i involved. ShaJfcin 
et. al. (Shn 67b) have solved this proolem by approximating 

with A. With an intelligent choice of A , e(A ) 

J-iir AiiA 

can be made close to and the difference V(§- can be 

rj! Sjjjv e 

evaluated by perturbation theory in I-IP basis. These corrections 
are referred to as spectral corrections and are expected to add 
roughlyc.-^lMeV per particle to the binding energies (Shn 67b) . 
Thus, binding energy' results would further improve when these 
corrections are talien into account. 
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It is expected that the results v/ould inprove if more 
terms are included in the expansion (3.20) of HP wave functions. 
Howrever , the situation is somev;hat ambiguous due to the following 
reason. It has been pointed out by Pal and Stamp (Pal 67) that 
there is a possibility of double counting of the contribution 
of high-energy states when HP wave functions are expanded in 
terms of oscillator states with large number of nodes and 
singular interactions are used. The strongly repulsive short- 
range part of these interactions gives rise to excitations to 
high-lying intermediate states. The contributions of these 
states are included in the calculation' of reaction matrix. How, 
when the HP wave functions are expanded in terms of harmonic 
oscillator states of n = 1, 2, 3 and still higher major shells, 
which are quite high-lying in energy, there is a possibility of 
double counting. Thus, 'one is left with two alternatives. First 
is, if one wishes to avoid the double counting only low-lying 
oscillator states should be included in the expansion of HP 
orbitals and one has to be satisfied with almost pirce oscillator 
states. This would not determine the tail of the wave function 
correctly. On the other hand, if more terms are used in the 
expansion (3.20), the HP wave functions arc determined more 
accurately but there is always a danger of double counting. 

Since inclusion of oscillator states of n = 1, 2 and 3 major 
shells appears to be a reasonable compromise- between these two 
alternatives, we expect that such a double counting would not 
cause any serious error in our calculations with effective Yale 


interaction. 
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In the oalculation of Sussex matrix elements (Elt 68a), 
as pointed out by ilavromatis (Mas 69) , the higher— order terms in 
the auxiliary potential are left out (see expansion 2.56)» a.he 

most important of these are the second-order Born term 7^7 and 

e 

the dispersion term ~ calculation of reaction 

^ ®n ®A ^ 

matrix elements mth Hsimada— Johnston (HJ) potential Kuo and Brown 
(Kuo 66) have observed that exclusion principle plays an important 
role in reducing the size of the Born tern for singlet states. 
However, they have also observed that there is a large contrib- 
ution to this term from the tensor part of the force for triplet 
states. The studies of Dahlblom (Dari 64) and Levinger et. al. 

(Ler 60) also show that tensor forces give large contributions 
to binding energj^ when calculated to second-order in perturb- 
ation theory: tlieir contribution being typically ten times larger : 
than the corresponding contribution for tlie pure central forces. 
Thus, there would be considerable improvement in the binding 
energy results v/lien this term is accounted for. Dispersion 
corrections are also expected to contribute significantly. 

In Tables (17.30 - 17.33) ^^0 and results of ; 

present calculation are compared with those of various other 
spherical HI and Brueohner-Hartree-Foclv (BHF) calculations. 
Different types of nucleon-nucleon interactions used in these 
calculations are also specified. It is seen that Yale and HJ 
potentials give very similar res-ults when the second-order term { 
in 7P is not talcen into account in the former. This is hardly f 
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Table IV. 32 

Compariscn of our results for the HI? proton single particle energies of 
with other oaloulations. All numbers are in MeV. 
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gable IV. 33 

Goinpari&..©n of our results for the H? single particle energies in ^^Gaon with 
other calculations. Numbers inside rounded brackets are the results obtained 
when second-order term in pseudopotential was included in effective Yale inter- 
action. Proton energies are marked with f, £11 numbers are in MeV. 
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surprising oecause of the similar nature cf these two potentials 
and also for tlie fact that they give same quantitative fit to 
the two-body data. Besides, in the calcula-tions referred in 
the table, reaction matrix elements from both these potentials 
were obtained using similar techniques. Por example, effective 
interaction from Yale potential v/as derived using koszkowski- 
Scott separation technique and an attractive pseudopotential 
¥/as added for those states where the potential is repulsive. 
Reaction matrix elements from HJ potential were also calculated 
using separation technique and the reference-spectrum method was 
used for the repulsive part of the potential. It is of some 
interest to compare the Sussex interaction results with those 
of potentials obtained from a fit to the phase shift data such 
as HJ, Yale^and labakin. It is found that Sussex interaction 
results for the binding energy are better than any other real- 
istic interaction results at least in first-order. Strikingly 
close similarit-y can be seen in the Sussex and effective Yale 
results when the second-order term in VT is included in latter. 
Very few spherical HP and BHP calculations have been carried out 
for and ^He. In the spherical HP calculations of Davies 

et. al. (Das 66) using a very simple velocity-dependent potent- 
ial the calculated B,E./A is -2.40 lieV, which is larger than 
the values obtained with Yale (-1.42 ^'eV) and Sussex (-1.99 MeV) 
interactions (see lables IV. 1 and IV. 18). Recently, Krieger 
(Krr-Preprint) lias reported HP calculations in a basis of 
cartesian harmonic oscillator wave functions using a velocity- 
dependent interaction especially derived for use in HP 
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calculations. For ^He, he obtains, B.S. /A = -0.825 MeV compared 
to our calculated values of -3*77 BaeV and -3.60 Me¥ v/ith effect- 
ive Yale and Sussex interactions respectively. 2hus, it appears 
that realistic interactions which fit the two-body data give more 

^ A. 

binding in He tnan the effective interaction of Krieger which 
fits the nuclear matter data. 

We have found that the binding of the lowest 
state increases very rapidly as the nucleon number changes* A 
local potential with reasonable strength vrould not give such 
deeply bound states and thus, this phenomenon is attributed 

to the non-locali'uy of the H? potential. Ihere is some experim- 
ental evidence from (e, e*y) experiments (Ami 64) that the 
binding energy of this level really changes from 35 to 6'C MeY 
as one goes from A = 12 to 28. However, the interpretation of 
the data is not unambiguous in view of the very wide peak cor- 
responding to Ko^re experimental evidence to settle 

this point conclusively seems to be necessary. 

HF is the first term in Brueclmer-G-oldstone perturb- 
ation expansibh end Second and higher-order corrections to it 
will arise from 2p-2h (t?/o partlclehtwo hole), 4p-4h, and higher 
excitations. Ihe large gap between the occupied and the un- 
occupied levels indicates. that these corrections would be small 
so that Hf solutions would be ■©lose to the actual solutions. 
However, in order to obtain good wave functions and also 
- improved -binding ener^, one must introduce 2p-2h type 
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correlations in tlie framework of the theory. The error 
introduced in choosing the harmonic oscillator wave functions 
as a "basis set can be considerably reduced by improving the 
wave functions at large distances to give proper tails. The 
results of our analyses also indicate that our approximation 
of assuming an average value of the oscillator parameter b (same 
for all states) nay not be good for malzing detailed calculations 
with Hi* wave fimctions. Ho?/ever, the uasic characteristics 
indicate that care must be excercised while treating the 
particle-hole 'tj-'pes of problems. 

* 5 Results for Sing-o Oscillator 
Configuratioii Approximation 

In this approximation, we have calculated the binding 
energies of ^He, ^^0 and ^^Oa. In the case when HF wave 

functions are replaced by pure oscillator wave functions, the 

4 - 16 . 

binding energies of doubly closed LS-shell nuclei (^He, 0, 

^^Oa.) are calculated using expression (5.73)- Wei^t-factors 
for ^^0 and are given in Tables (l¥.54 - IV. 55). has 

a partially filled p-shell and we have iised the following 
approximation to calculate its potential energy . 

P.E.(^^G) = ¥(ls. Is) + I [¥(ls. Ip) + W(lp, Ip)] , 

(4.6) 

where "¥(«, p) gives the contribution to potential energy arising 
from interaction of nucleons in the states a and p. Binding 
energy results for these nuclei are presented in Table IV. 36 
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Ta61e IJ/.34 
16 

Weight-f aotor.s for Og. We liave used the notation 
2S+1 for relative states. 

J 


State 

TJeight-f aotor 

^^S^(OOOO) 

3 

^^Sq(OOOO) 

3 

^^S3_(0001) 

9 

^^Sq(OOOI) 

9 

^^S2_(0010)l 

i 

^^S3_(0002) 

V 9 

1 

t 

/ 

^^Sq(OOIO)^ 

1 

i ^ 

5^3^(0002)1 

J 

L 

i 

i 

^■‘'P^(OIOL) 

6 

55p^(O10L) 

6 

55p^(010L) 

18 

55p2(010L) 

30 

^5d^(0200) 

1.5 

^5^^(0200) 

2.5 

^533^(0200) 

3.5 

5^0^(0200) 

7.5 

^5 3^(1000) 

1,5 

51sg(1000) 

1.5 


■Jtk A 
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TableIV.35 


Weigiit-f actors for 

Weight-factors for 


as in Tahle IV. 34. 


Fotation used is same 

states 


are same as those for 




states, 


State Weight-factor 


^^s^(oooo) 

3 

^^3^(0001) 

9 

^^3^(0002) 

15 

^^S^(OOIO) 

3.0 

^^S^(OOll) 

6.75 

^^s^(0003) 

15.75 

^^S3_(0020) 

1.125 

^^S^(0012) 

5.62 

^^S3_(0004) 

10.13 

^^S3_(1000) 

3.0 

^^3^(1001) 

5.25 

^^3^(1010) 

1.75 

^5s3_(1002) 

3.13 

^^3^(2000) 

11.12 

pp^(oioo) 

3 

53p^(0100) 

3 

^^P^(OIOO) 

9 

^5p2(0100) 

15 

^^P^(OIOI) 

9 

55Pq(0101) 

9 


State Weight-factor 


^^P^COlOl) 

27 

^^pgCoioi) 

45 

^^P^(0102) 

. 7.75 

^^Pq(0102) 

7.75 

^^P^(0102) 

26.25 

^^P2(0102) 

43.75 

^^P^(OllO) 

1.75 

^^Pq(OIIO) 

1.75 

^^^^(0110) 

5.25 

55p 2 (0110) 

8.75 

^^P-, (0111) 

2.25 

55p^(0111) 

2.25 

^5p^(011l) 

6.75 

55p2(0111) 

11,25 

^^P^(0103) 

5.25 

^5?j^^C103) 

5.25 

55p^(0103) 

10.75 

^^p^Coios) 

26 . 25 

^^P^(llOO) 

2.25 

55Pq(1100) 

12.25 


Gontd 


gable IY.35 (Gontd.) 


C? - 1 - „ -i- 

D b a ti 0 

7J e i gli t - f a c 1 0 r 

Q J- * 1 - ^ 

0 -j uzi 

77eigiit-f actor 

^^P-, (HOC) 

6.75 

13 , 

13(0210) 

1.46 

^^P^dlOG) 

11.25 

^^1^0210) 

3.13 


2.25 

^5?.. (1200) 

a. 

1.12 

53e (1101) 

2.25 

-5? 2(1200) 

1.38 

55pdllCl) 

6.75 

1%5(1200) 

2.62 

55? 2(13-C’1) 

11.25 

5112(1200) 

5.62 

^^I>^(0200) 

3 

^•"--(0300) 

5.25 

^0200) 

5 

53 

I'2(0300) 

11 .25 

15 , 

1)3(0200) 

7 

^^13(0300) 

1 i7 7 C 

51 , 


33 . 


DgCOlOO) 

15 

j:-£ 0 500 j 

20.25 

^^]D^(0201) 

5.25 

^^^3(0301) 

!? • 2 3 


8,75 

^^2(0301) 

[] — L .>23 

^^13(0201;' 

12.24 

^^3(0301) 

15.75 

5^D„(0201) 

26 .24 

(0301) 

20.25 

13 / 

D^(0 20 2) 

13 / 

K- (020 2; 

3.CS 

-j 3J, H" 

^0„(O4OQ) 

2.63 

5.13 

^ •rr ^ 

~%( C 400 ) 

3.32 

1513(0202) 

7.20 

^5?^ (0400) 

^.13 

51d2(0 2C2) 

15-42 

% A 0400) 

10.13 

^^2^(0210) 

0-62 



) ^^I!2('52iC) 

1.04 






gable IV. 36 

Binding energy per particle of Tarious nuclei 
obtained in tbe single oscillator 
configuration approximation. 


Interaction 


Sus 

sex 


Yale 

nucleus 

^Heg 


^8 

'^Ooa 


b(l) 

1.4 

1.8 

1.8 

2.0 

1.76 

B.E./A (MeY) 

0.33 

-3.72 

1 

* 

0 

1 

-3.14 
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for tlie values of 7) which give r.m.s. radius close to experim- 
ental values. Eosults obtained with effective Yale interaction 
are also given in this table. In the binding energy calculatior 
of ^^Oa with Sussex interaction and for b = 2. OF we find that 
contributions to potential energy of S, Pj D and P states are 
93-42' j 0,L':>2 5 5-52 and 0.552 respectively. The G— state 
contributes alnost' negligibly. P-state raatrix elements are 
large but they tend to cancel each others contribution to 
potential energy . On the other hand, although P-state relative 
matrix elements are relatively smaller, they occur with large 
weight-factors and make a net contribution to potential energy 
comparable to that of P states. Thus, in any precise calcul- 
ations P-state should be included while G-state can be neglected 

16 

Single particle neutron and proton energies of 0 

obtained with Sussex and Yale interaction are given in Table 

IV. 37 (see Section III. 9 for formulae used). It is seen that 

Sussex levels are more spread out for smaller values of b. Single 

40 

particle neutron energies for Ca are given in Table IV. 38 alor.^ 
with the experimental values. On comparing these energies with 
self-consistent 11? single particles energies (see Table IV.23)» 
it is seen that the effect of self-consistency on the single 
particle levels is to increase the binding of the occupied 
states. Same conclusion can be drawn on comparing the ~ 0 leve? 
given in Table IV. 37 with tliose given earlier (see Tables IV. 1 
and IV.19) . 
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Cohen et. al. (Con 63a, 63b) 
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Table IY.38 


40 '' 

Neutron single particle energy levels of ^ GagQ 

for b = 2.0 


singl e 

particle 

state 

SUSSEX 

(IleY) 

Experimental®' 

Id^ 

p 

- 21.8 

-21,9 

2 s^ 

2 

- 21 . C 

-18.2 

Id, 

p 

2 

-IS. 85 

- 15.6 

If^ 

2 

- 5.06 

- 8.4 

2 P 3 

- 6,C3 

- 6.3 

Ifr 

5 

2 

- C.99 

- 4.3 

2P]_ 

- 3.07 

- 2.9 


a. Elton et. al. (Sin 67) 
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On comparing the single oscillator configuration 
approximation results with those obtained in El calculations, 
we find that the approximation of using pure oscillator wave 
fxinctions is not too had as far as the estimation of binding 
energies and single particle energies are concerned. 



GHAP2BR V 


COIRIGURATION MIXING CiliCULA'IIONS f/IIE 
HARTRBE-FOCK WAVE PDlfGHONS 

V.l Details of tlie Galoulations 

In this Chapter we shall describe the restilts of 
configuration aizing calculations, with Hartree-Eocli; (HP) wave 
fxmctions^ carried out in the spirit of a standard shell model 
spectroscopy calculation for the nuclei and Ihe 

main assumption of the sheU. model is that the nucleons move 
independently of each other in an average potential well, com- 
monly taken to be of harmonic oscillator or Woods-Saxbn type# 

She more fundamental way of obtaining the average potential 
field of the shell model, starting from free nucleon-nucleon 
interaction, is the HP method. We have seen in Chapter IV that 
the HP orbitals for the occupied levels of low-mass nuclei can 
be well approximated by pure oscillator wave functions but the 
two differ appreciably for the unoccupied levels. Since the HP 
orbitals are obtained in a more fundamental way, one would be 
tempted to believe that they should be closer to the true single 
nucleon wave functions compared to the pure oscillator wave 
functions. She single particle wave functions and energies 
used in our' calculations have been obtained by solving the HP 
equations seif-consistently as described in Chapters III and I¥ 
Each unloiown EF sin^e particle orbital is expanded in terms 
of pure oscillator Y/ave functions 


( 5 . 1 ) 
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Here a and ji specify the angular and isospin quantum numbers of 
the Hartree-FocI: (HP) and harmonic oscillator (HO) single 
particle states; a = (al/omr), p = (nijiiiT) , where a and n are 
radial quantum numbers, The HP field is constrained to be 
spherically symmetric (see Chapter III)j and also the forces are 
assumed to be charge independent so tliat the summation in equ- 
ation (5.1) is restricted only to the radial quantum number n. 
Three oscillator states of b- =1, 2 and 3 major shells are used 
in each expansion. Thus, equation (5,1) becomes, 

n=l 


7atji 


0?(t3) f. 


nvimx 


(5.2) 


In a spherical field the angular partsof the single particle 

m 

wave functions are the vector spherical harmonics j.defined 

in equation (2.14) and thus, more explicitly, the wave functions 
and 9 can be written as 

T-'Q Dl ^ T 


^nijmr = 


(5.3) 


where E and% specify the radial (see expression 3*19) and 
isospin parts of the wave functions respectively. Using 
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expressions (5.3) > ec^uation (5.2) can be rewritten as, 


Hi* 


S 


Ctj = S 


n=l 


(5.4) 


The low-lying excited states of and are 

obtained under tlie assumption that the last two nucleons move 
in a HI field and interact with each other throu^ some 
residual interaction. Ihe remaining (A-2) nucleons take part 
only in producing the average field, She energy levels are 
obtained by constructing and diagonalizing a Hamiltonian matrix, 
including all the degenerate configurations having the lowest 
possible energies. file matrix elements are <( a^JT” » 

where a, p, y, 6 specify the HP single particle states. These 
matrix elements are evaluated analytically using expansion (5.2), 

p/ 






(5.5) 


Writing explicitly. 


<C^tl^l» ^2^2^2* 


. — ■ — , 21 -| * / 23 .#^ ^ s ^ i i 

^ ^ ^a^^^l^l^ ^ap^^2^2^ ^^^^1^1^ ^a^P^^2^2^ ^ 


/ / 1 


n^ t 
^2 


X ^^l^^l* ^2^2^2‘ ^2^^2V^'^>H0 ' (5*6) 
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Thus, each of the El matrix element is expressed as a sum of 
several two-body matrix elements in harmonic oscillator basis 
with proper weight-factors (products of different G-coefficients) 
There are eighty-one terms to he summed in eq_uation (5.6) when 
3 states are used in expansion (5-2). The expansion coef- 
ficients are obtained self-consistently and the results have 
already been presented in Chapter IV. The two-body matrix 
elements of the interaction in harmonic oscillator basis have 
been calculated using equation (3.26). A maximum of relative 
h = 2 states has been used in our calculations. Single part- 

icle energies are added to the diagonal terms. ¥e have carried 
out the calculations with both the IH single particle energies 
calculated in Chapter IV» and the experimental single particle 
energies. The l£.tter approach has been followed by most of the 
authors (Kuo 66; Clt 68). It has got the advantage that it 
takes into account the interaction of the last two neutrons 
with the closed core to all orders, but the disadvantage is 
that it may obscure the shortcomings of our calculations which 
might be obvious otherwise. 

The following two nucleon-nucleon interactions have 
been used in our calculation: 

(i) Effective Yale interaction (see Section II. 2) 

(ii) Sussex interaction (see Section II. 3) 

The prescription used to obtain the effective Yale interaction 
from the Yale potential (Shn 67a), and the matrix elements of 
the Sussex interaction (Bit 68a) in an oscillator basis 



directly from experimental phase shift data ha's been described 
in Chapter II . 

In an ideal case, the shell model calculations should 
be carried out in an infinite dimensional- configuration space 
when the free nucleon-nucleon interaction is used as residual 
interaction. Iliis is, however, not practicable and the shell 
model calculations can only De carried out in a limited config- 
uration space. Tile neglected shells renormalize the interaction 
in the space of the shells considered in one's calculation. Ihe 
problem is then to find out an effective interaction which would 
give the same results, when used in a limited configuration 
space, as the true interaction would give in an infinite dimens- 
ional configuration space. Several authors have worked in this 
direction of establishing a relationship between free nucleon- 
nucleon interaction and effective interaction. We mention here 
the work of Kuo and Brown (Kuo 66, 67) dealing with Hamada- 

Johnston (HJ) potentialj and of Clement and Baranger (Clt 68) 

18 

dealing v/ith labakin potential. In their calculations of P 

1 R 

and 0, these authors find that there are two important cor- 
rections arising from neglected configurations to the free 
nucleon-nucleon interaction. Pirst is the second Bom correction 
in the matrix elements of triplet even tensor force. Ihese are 
the corrections arising from the excitation of two valence 
nucleons to intermediate states of very high-energy. Ihese 
corrections mainly affect the 'f = 0 matrix elements- In deriv- 
ing the effective Yale interaction from the Yale potential, 



Shakin et. al. (olm 67a) h.a.ve oaJ-culated the second-or-ler Bom 
corrections "by combining Moszkowshi-Scott separation teclinique 
and kuo and Brom method (Kuo 66) of taking intermediate states 
as plane waves amd angle averaged Pauli operator. No such 
corrections are included in Sussex interaction. 

'Ihe second important correction arises due' to the 
excitation of the panticles from 0 core to the intermediate 
states, which lie close in energy to those included explicitly 
in the diagonaliza,tion. Such corrections are known as core- 
polarization effects. Kuo and Brown, and Clement and Baranger 
have calculated the admixture of 3p-lh (three particle-one hole) 
states to the usual 2p state by using second-order perturbation 
theory in harmonic oscillator basis. Ihey find that these cor- 
rections, although relatively smaller than the Born cor- 
rection, are quite important for^ = 1 states particularly for 
1 = 0"^ states. ^=0 states are, however, almost unaffected. 

Ko such corrections have been included in our calculations for 
either of the two interactions used. 


The configuration space used in our calcurations for 
A ■= 18 nuclei consists of (id^^^ shell model states. 

The ground state configuration for both these nuclei is 

X8. 

(id^^g)^- The level spectrum that one would observe in 0 is, 


(ld^y/2) • 

J = 0"^ 


1=0"^ 

(Idj/g aaiL/g) : 

J « 2 ^ 


, 2 ^ 4 -^ 
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111 these states have total isospin 't' = i, ihe following‘X'= 0 
levels (besides the T* = 1 levels mentioned abov^ would, be 
observed in the case of 

(ld5/2)^ : j = 1+, 3+, 5+ 

(2s^^2)^ ; J = 1+ 

(11^/2 : J = 2"^, 3^ - 

fhe ldj ^2 state, v/hich lies at 5MeV in ^"^0, has not been 

included in our calculations mainly because of computational 

difficxilties, She HI* single particle parameters used in our 

IS 

calculations of 0 should have been obtained from self- 

17 

consistent HI* calculations for. '0. Ihe observed grotmd state 
17 

spin of '0 is J = ■“ . The HF method described in Chapter III 

17 

would only describe an intrinsic state of 0 and the ground 
state with spin ^ can be obtained by projecting out good 

angular momentum states, lo avoid such complexities, we decided 

to take single particle parameters from sclf-consistent HP 

In 1 R 

calculations for 0 or 0. She spherical HP calculations 

give the J = O’*" ground state of these nuclei. In taking the 

1 

single particle parameters from 0 HP calculations the contrib- 
utiore of the interaction of outer two neutrons with the closed 

core to the self-consistent average field and wave functions 

IB 

are missed. In 0 HP calculations, part of the interaction 
between the last two neutrons is already included in the calcul- 
ation of self-consistent field and thus, the residual inter- 
action left over is not just the free nucleon-nucleon inter- 
action but gets modified. However, wo think that these 



corrections would oe small and average self-consistent field 
would not be very different for ^^0, ^"^0 and . In all our 
calculations tlie last two neutrons in ^^C*, and one neutron and 
one proton in liave been treated on eq’aal footing. It is 

assumed that the self-consistent field is the same for neutron 
and proton so thei-i; same wave functions (neutron wave fimctions) 
describe the motion of both the particles outside the closed 
core. This is, iiov;ever, not a very serious restriction as we 
have seen in Chapter IV that coulomb force has little effects 
on the single particle HP wave fractions of light isaclei. 


10 . 


The configuration space used for the calculations of 
Be levels consists of Ip^ /o and Ip-, shell model states. 


V2 ^^1/2 

The ground state configuration for this nucleus is (Ip^^^) . 

All low-lying levels of ^^Be woiild have 'T' = 1 and the following 


spectrum would be observed} 


(lP3/2)^ 

j 1=0"^ 

(IP3/2 ^^1/2^ 

11 

H 


: 1 = 0"^ 


V.2 Results and Discussion 

Results of our calculations obtained with effective 

18 

Tale interaction are shown in Pigures (V.l - V.2) for 0 and 
^®P respectively. In these figures, YA and YB refer to calcul- 
ations done with single particle wave functiona obtained from 

1 fi 

HP calculations for 0 without coulomb and centre-of-mass 




YA YB YC YD YE YF EXPT. 

F!G.Y-2 Positive parity states of obtained witb effective 
Tale interaction and Hf wave functions. Levels, witn 
T *a 1 are specified wltli braol^eted spin* All otber 
levels have T » 0. 
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corrections, and using HI* and experimental single particle 
energies respectively (see Tables IF. 11(b) and IF. 5). Cases IC 
and YD refer to calculations done using uncorrected (without 
coulomb and centre-of-mass corrections) IIP wave functions for 
and with hf and experimental single particle energie® respectively 
(see Tables IF, 13(a) and IF. 5). Cases YE .and YE refer to calcul- 
ations done using experimental single particle energieSt and wave 
functions obtained from HE calculations ?/ith coulomb and centre- 
of-mass corrections for the nuclei ^^0 and respectively 
(see Tables IF. 11(a) and IF, 13(a)). All calculations with 
effective Yale interaction have been carried out for b = 2.09E 
since for this value of b (=Yb/m^l>) we could obtain minimxun 
total energy for 0 in self-consistent HE calculations (see 
Figure IF.l) . For consistency, calculations for Sussex inter- 
action have been done for b = 2. OF. The experimental single 


particle energy difference between and 2 s^y 2 state was 

taken to be 0.87 licF, as obtained from experimental levels of 
^'^O. This is not justified for ^^F for which an average of the 


neutron and proton single particle energies for 2 s ^^2 states 
should be talcen. However, the discrepaney would be small. The 
experimental levels for and ^®F, shorn in Figures F.l and 
F.2, have been tal'-eh from references (Hea 64) and (Poi 65) 
respectively- The absolute ground state energies are determined 
as follows. Using the binding energy data of Koenig et. al. 


(Kog 62), the ground state energies of valence nucleons are 


taken to be 
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- B.E.(^S) + 2 B.S.(^'^O) = -3.90MeV for 

and 

-B.S.C^^E) - B.E,(^S) + B.B.(^'^O) + 3.S.(^'^P) = -5.00 lie? 

for 

As seen fron the figures, except for the 4"^ ('X=l^ states 
the calculated levels are in very poor agreement with experim- 
ental results in all the cases. In particular, the ohserved 
ground state of is a 1 "^ state wherea,s in all our calcul- 
ations we get 5"^ state as lov/est state and first 1 "^ state lies 
much higher. However, the fact that very similar results are 
obtained in various cases supports our assumption that HP field, 
if constrained to remain spherically symmetric, would not he 
very different for ^*^0, ^"^0 and Results obtained with 

effective Yale interaction using pure oscillator wave functions 
and experimental single particle energies are shown in figures 
(¥.3 - ¥.4) forC =1 and'C'^0 states respectively. In these 
figures, HOI and 102 refer to calculations done with harmonic 
oscillator (HO) wave functions in configuration space consisting 
of (l<i 5/2 ^^*^ 5/2 ^®l /2 ^S/ 2 ^ states respectively. 

Ihe importance of including the 1 ^ 2/2 is obvious 

from the figures. It lowers theT=l levels, thereby improving 
the agreement v;ith experiments- and also the first 1 level in 
^®f is pushed below the lowest S"*" state. In figures (¥.3 - ¥.4), 
we also compare the calGulatedC=l and"0 levels obtained with 
Hf wave functions in case YB with the results of Kuo and Brown 
(referred to as KB), and Clement and Baranger (referred to as 



ENERGY CMeV ) 





ENERGY CMeV) 



FIG -2-4 Comparison of calculated- ilf and HO positive paritjr 

lovels of for effective Yale interaction idtH the 
nMiats of Kuo and Brotm (KB), and Clement and Baranger 
(0^2). Levels with Bracketed spin have T * !• All 

oulex levels have T » 0,* ■ 



GB2 when secoiid-ordex Born tern is included) . Both these calcul- 
ations were done in a configuration space consisting of 
^^ 1/2 ^^ 3 / 2 ^ hsirnonic oscillator states and using experimental 
single particle energies. E3 and GB2 levels shown here include 
the oontrihution of second-order Bom term hut not the core- 
polarization effects. On comparing the HB spectrum obtained in 
case YB and the ilOl spectrum (see Figures V.3 - Y.4) we find 
that other parameters being held fixed, IIF and HO wave functions 
do not give verj different results. 

In Tables ?.l and Y.2, v/e present the calculated two- 
body HO and HP matrix elements in the s-d shell obtained in 
various cases. Results of several other calculations are also 
given for comparison , Federman (Fen 55) and Arima (Ara 68) 

*^=1 matrix elements have been obtained from best fit to experim- 

13 rr' 

ental levels of 0. Arima T=0 matrix elements are obtained 

1 Q IQ IQ 

from best fit to experimental levels of 0, F and He. The 
values obtained in two different cases, first when all but the 
1.7MeY, I"*" state in ^^Fwere included from the fit and second 
when this state was also included are referred to as A1 and A2 
in Table Y.2. It can be seen that no impressive agreement of 
the calculated ilF matrix elements exists Y/ith any one of these 
calculations. In particular, the J=0'^ HP matrix elements are 
considerably smaller than the results of other calciilations , 
However, Federman-Talmi (PT) and Arina have obtained positive 
values for J=3'^ and i"*” (r=l) matrix elements. In oxir calculations, 
j=5'^ ^'^=1] HF matrix element is always positive although much 
smaller in magnitude. 
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Tile results obtained with. Sussex interaction are shown 
in Figures (7.5 - V.6) for ^^0 and respectively along vdth 
the experimental levels. In these figures, SA and SB refer to 
calculations done using uncorrected HP v/ave functions for ^*^0 
(see Table I¥. 27 (b) ),, said HP single particle (see Table IY.24) 
and experimental single particle energies respectively. HOI and 
H02 are the results obtained v/ith pure oscillator wave functions 
in a smaller (ld^y'2 and larger (ld^^2 ^^iy^2 ^*^3/2^ 

configuration space respectively. Glement and Baranger results 
obtained ¥/ith bare Tabahin interaction without second-order 
corrections (referred to as CBl) are also given for comparison. 
In Table ¥.3> calculated two-body HP matrix elements for Sussex 
interaction in tie s-d shell have been compared with Clement 
and Baranger results for Tabakin intersection (031). 


10 

The calculated energy levels of Be for effective Yale 
interaction are shown in Figure ¥.7 along with the experimental 
data. In this figure, XA and XB refer to calculations done 
using HP wave functions for (see Table I¥.12), and HP single 

particle (see Table I¥.5) and experimental single particle 
energies respectively. XG is the result of calculation using 
HP wave functions of ®Be (see Table I¥.9) and experimental single 
particle energies* Again, HOI and H02 are the results of calcul- 
ations using oscillator wave functions with experimental single 


particle energies in configuration space consisting of (IP3/2 
IPl/g) and (IP3/2 ^^ 1/2 ^^ 5 / 2 ^ states respectively. The experim- 
ental ener^ separation between lp^^2 ^^1/2 


ENERGY 


> 



FiG-ST-S Comparison of calculated HP and HO positive parity ■ 

18 a 

levels of 0 for Sussex interaction with the results 
of Clement and Baranger (results without second-order 
Bom correction are referred to as CBl) , and expericea 
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FIG -2 -6 Compariion of oaloulated Bf aiMi 10 pooittve parity 
( r w 0) l®T«ls of f«je Su»a«3C fatesaotion witl 
tlo reatilti of Olaaaat amSL Baxaagox (CBl) , and 
exporimtfata. ’ ' . 
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to be 6MeV. All calculations have been done for 13=2.09?, 7e 
find that the csilGulated difference between the ground state 
and first 2"^ state is much less compared to experimental 
separation. Eo^vever, the calculated energy separation between 
the first and the second 2'^ states is 5.85 Met in case IB 
compared to experimental separation of 5.9 Met. 

In all our calculations we find that calculated H? 
spectra are compressed compared to experimental levels. She 
poor agreement vri.th the experimental results may be due to the 
nature of the nucleon-nucleon interactions used in our calcul- 

I 

ations. As seen in Chapter IV, both these interactions give 
3--4MeV less binding when used in spherical HP calculations. 
However, spin-orbit splittings in vanious nuclei have been 
satisfactorily reproduced with both these interactions (Shn 67b: 
Elt 67). Thus, it appears that the spin-orbit part of these 
interactions is reasonable but the central part is weak. In 
yet another entirely different approach, Gimye has calculated 
(Gue 68b, 69) the low-lying states of various odd and even mass 
light nuclei by projecting out good angular momentina states 
from an intrinsic deformed HP state. He has carried out the 
calculations for 2s-ld shell nuclei (Gue 68b) with effective 
Yale interaction^ and for Ip shell nuclei (Gue 69) with both 
effective Yale and Sussex interactions. In all the eases the 
projected HI spectra for even mass nuclei are found to be compr- 
essed compared to the experim^tal levels. For example, in the 
case of the calculated separation between the ground state 



and first 2 ^ state is 0.51 MeV compared to the experimental 
value of 1.93 ivIeV, llie effect of increasing state matrix 
elements arhitrarily oy over 30 percent so as to reproduce the 
experimental binding energies of various nuclei is to spread 
out the projected spectra. Since the results obtained by 
projection method a.re close to shell model results (Eeh 58), 
the compressed spectrum obtained in the present shell model 
calculations in ItlF' basis seems to be consistent with the find- 
ings of Gunye. 

It should be emphasized tha.t the poor agreement of our 
results with the e^cperimental data does not imply the shortcoming 
ings in the HP bc^sis since we have not made any attempt to 
calculate the important second-order corrections in both the 
interactions arising from core excitations (besides the Born 
term for Sussex interaction), To the best of our jfcnowledge, no 
attempt till now has been made to use self-consistent single 
particle energies and wave functions obtained with realistic 
potentials in the shell model calculations. Kuo and Brown 
(Kuo 66) have calculated single particle energies in '0 with 
their realistic potential and obtained qualitative agreement 
with observed energies, but they did not use calculated single 
particle energise in obtaining the spectra of A=18 nuclei- 
Our calculation is a first step towards a more complete shell 
model calculation in the HP basis. In fact, it has been shown 
tlie.t low-lying states of are abundant in deformed components 
(Bron 66: End 65) and the assumption of two nucleons outside 
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16 

a closed 0 coxo is not very good. 2iie transition rate froa 
second 0 to first 2"^ state is also very large and suggests 
that these stares should be mainly composed of deformed comp- 
onents. However , Kuo and Brown (Kuo 66) have satisfactorily 
reproduced the enorgy levels of these states by using simple 
two particle model mth p-h corrections. It is expected that 
such core excitation corrections Y/ill spread apart the calcul- 
atedT=l HI spectrum by depressing the lower states and raising 
the higher levels and thus, the agreement with the experiment 
would improve. Studies of Kuo and Brown, and Clement and 
Baranger show that in the harmonic oscillator basis most imp- 
ortant of these corrections are 3p-lh- corrections arising from 
excitation energies of 2ft£U. Contributioiis of higher levels, 
which involve excitation energies >2h.u>, a.re relatively unimp- 
ortant. However, this may not be true in I-H basis. Bassichis 
et. al, (Bas 69) have made explicit calculations for the second- 
order energies and probabilities of 2p-2h components in the 
ground state of ^he and ^^0 using' perturbation theory in HI 
basis. Ihey find that the relative phases (signs of expansion 
coefficients) of the HI wave functions of unoccupied states 
play an important role in deciding the importance of a part- 
icular level. It is quite likely that excitations to some of 
the high-lying levels may be more important than to low-lying 
levels contrary to the predictions of simple shell-model 
picture in harmonic oscillator representation. However, no 
conclusive remarhs can be made about the magnitude of such 



corrections in HP basia unless explicit calculations are 
made , 


In another calculation carried out for the 
residual interaction 'between the last two neutrons in 
was talien to he of the form (a + ^'^•'02) e^cpC-rV^o), Iwo- 
hody matrix elements in the s-d shell were calculated in a 
basis of HP wave fimctions and we tried to determine the 
parameters of the potential v/hich would give a good fit to 
the experimental energy levels of But no set of 

parameters could satisfactorily reproduce the experimental 
levels of -^ 0 . Thus, it appea.rs that a simple residual 
interaction of the form described above is not suitable 


for a HP basis. 



iUPPSlTDIX 


Tlie local potential; used in our calculations of 
single particle energies, is of the following form. 






so 


where 


(A.l) 


' ^ 1 ^ 

f(r) = 1 + exp(— V 


■I 


f 

/ SO'' • 

ezp(--- 


j 


r-R, 


-2 


) <;l + exp( — ) 


g(r) = — 


SO 


1 


a 


so 


r 


1 


li(r) = ^ 


_i.( 3 _ 

OT> V pj 

Rq 


for r^E^ 


for 


(A. 2) 


Y,^ = Central nuclear potential, 

Y = S’Din-orhit nuclear potential, 

SO 


(|j°cr) . = h 


for 3 = 


= - il+1) for 3 = 


(A.3) 


(xiix) 



th 


^'x radius and diffuseness paranexers of tils x 

part of the Potential and H = r 

X X 

<^ore mass and its charge » 

P'S are the particle reduced mass and its charge. 

The radial W3,ve functions s.re the solutions of the 

following Schrodinger equation, 


The "binding energies of the single particle levels are fed in 
and the strength parameters Y-^ and Y^^ are obtained by the 
method of successive iterations till the values converge. 

The v/ave functions are normalized to unity. 

A separable form of the local potential is given by. 


7(r)'4’(r) = H(! r^r^ i)v|; ( rO dr^ . 


(A. 5) 


H(j’r-'?^}) is given the Gaussian form, normalized to unity, 


H(ir^? }) = .a^exp) 1, 

i ‘^nl j 


(a. 6) 


where a^ is the non-locality range parameter. An approximate 
relationship between the wave function calculated in a 

potential of the form, given in equation A. 5, and the local 
wave function j is given by. 


■2 2 V ' o 

2Prtr ^ 0)g(r)-Z_^Zj_^e^hCr)] (a.?) 
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APPEIQIX 3 


It has osen mentioned before that the tv;o-body 
coulomb and ceiitre-of-nass matrix elements used in the HP 
calculations v/itli effective Yale interaction were calculated 
with a maximum relative 1 ( see Section lY.l) . This has 
caused some error in our calculations since higher I; values 
are also quite important for both of these matrix elements. 
The nevi coulomb and centre-of-mass matrix elements in the 
configuration space consisting of '^^3/2 ^1/2 

oscillator states of n = l, 2 and 3 ma-jor shells were 
calculated with maximum relative t =10 . The revised HP 
results obtained by using these new matrix elements are 
shown in Table 3.1 for the nuclei "^He, ^Be, ^^0 and ^^0- 
On comparing with the older results (see Table lY.l) we 
find that the effect of truncation at (j=2 on most of the 
static properties of closed shell nuclei is small. Hovrever, 
the wave functions are affected^ mainly in the tail region. 
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IPPMDIX 


n 

‘sj 


I'u is assumed that Sussex inter action V can he broken 
up into central, spin-orbit and tensor parts which are tensors 
of rank zero, one and tv/o respectively in the sca.lar product 
^ Here refers to the spin variables and 

to the orbital motion. (Thus, vie can write, 


10=0 


. 2 

I :=0 

= 7(0)+V(l)+V(2) (G.l) 



and the matrix elements of V(k) are 


<n' S ?j J I V( k) |nS '> = < n* S j| | nSi?J > 


= (- 1 ) 


J+S+ € 


s t J] 




i 


<s||s^’'^ II S>X 


S kl 


X <n' 0 n I >. (0.2) 


Smamiiig oveu Iz on botli the sides ¥/e obtain j 


<n‘ St'j ?|nS&J > =^(-1)-^+®^''' ] p i<S j| S> X 


US k 


X<Yt' ||v'=‘> |nt>. (0.3) 


(xvii) 



Multiplying Doth the sides hy 




V 


\l 


k 


and summii 


ig over J we get 




V J| 


S k 


S'i! J\Y\nshy 


■"N 


/ fs t j] [s 6' j 


5:1 WWt. 

kg - I S ic I i t 

J 


rK 


X <S II n s> <n' t' H ni.y 


= Z^Mc'<S 1! s('^)|l S> < 5 ..!' 11 r^'^Wniy 

k 

= <^sil ji niy . (G,4) 

Ihus, in equation (0.4) the \mcoupled matrix elements of any 
particular component are expressed as a sum of matrix elements 
of the total potential V in different channels. The coupled 
matrix elements of V(k) are obtained by substituting equation 
( 0 . 4 ) in equation (0.2). 


file off-diagonal matrix elements of the Sussex 
interaction V, with radial quantum numbers differing by 
more thaii one, have been calculated using the follo?;ing 
expression, given by Elliott et. al. (Elt 68a), 


(xviii) 




Kie relation (G.5) connects the matrix elements wi 
n^ -n = 2 tc those with n^ -n4 1 and is true for the matrix 


elements of any general interaction Y in an oscillator basis. 
It is assured that interaction Y is raomentum independent so 

p 

that comnutshor [r , V] vanishes. 
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